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Abstract 



In this article, I study the diffusive behavior for a quantum test particle interacting 
with a dilute background gas. The model I begin with is a reduced picture for the 
test particle dynamics given by a quantum linear Boltzmann equation in which the gas 
particle scattering is assumed to occur through a hard-sphere interaction. The state of 
the particle is represented by a density matrix that evolves according to a translation- 
covariant Lindblad equation. The main result is a proof that the particle's position 
distribution converges to a Gaussian under diffusive rescaling. 

(N 
>. 

ov 1 Introduction 

m 

CSl | Boltzmann's transport equation generates an idealized dissipative dynamics for the one-particle 
phase space density of a gas of interacting particles [21 H] . This dynamics is relevant for a dilute 
gas in which the essential information for the interaction between gas particles is reduced to 
two-particle scattering data. Boltzmann's equation effectively governs the convergence of the 
one-particle phase space density to the distribution corresponding to the thermal equilibrium 
state. The linear Boltzmann equation is similar in character but governs the dissipative evo- 
lution for the one-particle phase space density of a population of (non self-interacting) test 
particles interacting with a dilute background gas in equilibrium. By normalizing the density 
by the number of test particles, the linear Boltzmann equation serves as a Markovian master 
equation describing the noisy trajectory of a single test particle in the reservoir. A mathemati- 
cal derivation of Boltzmann's equation from a microscopic model for a gas of hard-spheres has 
appeared in the article [30J, which verifies a mathematical conjecture in [16] concerning what is 
referred to as the Boltzmann- Grad limit. A derivation of the linear Boltzmann equation from a 
microscopic model for a test particle interacting with a background gas follows from the same 
analysis and is discussed in the review [37] . Starting from a quantum version of a Lorentz gas, 
a linear Boltzmann equation has also been derived in a low density limit in [10J. 

The linear Boltzmann equation is a classical model treating the test particle as a point. A 
motivation for having a quantum dissipative dynamics rather than a classical one is to have a 
model in which the respective time scales of dissipation and decoherence for a quantum particle 
interacting with a gas are both apparent. The Lindblad equation, serving as a quantum linear 
Boltzmann equation, should be determined by the essential features of the gas (e.g. spatial 
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density, temperature) and the scattering behavior for the interaction between the test particle 
and a single gas particle. 

There have been various quantum master equations proposed in the physics literature to 
model the dissipative/decoherence influence of a background gas acting on a test particle [25], 
EES G3 EH1 [23] . With the microscopic picture for the test particle interacting with an ideal gas 
as inspiration, the derivations have all been heuristic or based on an analysis of the effect of 
a single scattering with a gas particle. The first heuristic derivations were [25] and [IT] where 
the authors considered limits in which the test particle becomes infinitely massive compared 
to a reservoir particle. The limiting dynamics in those cases are transient in nature, since 
by neglecting frictional effects, they predict linear growth in the mean kinetic energy. The 
first proposed model to include energy relaxation was in [9]. The model was later refined and 
given clearer motivation in [381 139] and subsequently [231 121] . A comprehensive review of the 
literature can be found in [ID]. 

In this article, I begin with the quantum linear Boltzmann equation discussed in [30] for the 
special case in which the interaction between the test particle and a single reservoir particle is a 
hard-sphere potential and the reservoir particles are distinguishable. The main result is a proof 
that the particle diffuses spatially in the limit of large times. My analysis yields a closed form 
for the diffusion constant in terms of the scattering amplitude for the hard-sphere potential. 
The techniques here are an extension of those in [3], which require the Lindblad generator and 
other operators involved be bounded. 

I will mention a few recent mathematical results regarding diffusion in quantum models 
before proceeding to define my model. The articles [HI [12] prove diffusion for a quantum 
particle in a random potential (quantum Lorentz gas) in a weak coupling limit. The article [26] 
shows diffusion for a quantum particle in a tight-binding approximation (living on Z d ) whose 
wave packet evolves according to a stochastic Schrodinger equation. The results of [26] are 
extended in [19] to handle a weaker form of stochasticity in which the random potential forcing 
the particle is periodic rather than evolving independently at each lattice site. Also for tight- 
binding models, [8] and [7J prove diffusion for microscopic models of a particle interacting with 
an environment. The environment in [8] is a lattice of non- interacting reservoirs that are taken 
to interact with the test particle in a weak coupling limit. The result in [7J treats a quantum 
particle interacting with a free Bosonic field in the limit that the particle is massive (weak 
coupling). 

In Sect. [21 I present the model, state the main theorem regarding the diffusive behavior of 
the model, and give some background information. Section [3] discusses the technical questions 
regarding the existence of the dynamics I consider. Sections 010 build up the technical tools 
for the perturbation argument that I employ. Finally, the proof of the result is contained in 
Sect. El 

2 Model, results, and discussion 
2.1 The model and the main theorem 

The state for the test particle at a given time t G IR+ will be described by a density matrix 
pt G £>i (L 2 (IR 3 )) , where £>i("H) denotes the space of trace class operators over a Hilbert space 
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%. The state p t evolves according to the Lindblad equation 

j t p t = C( Pt ) = -i [H, Pt ] + *(p t ) - ~{**(I), pj, po = p, (2.1) 

— * — * 

where the Hamiltonian H = H* is a function H{P) of the vector of momentum operators P 
with Pj = -ifigf: for j = 1, 2, 3, * is a completely positive map defined below, and 

is the evaluation of the identity for the adjoint map \E r *. The term — l[H, •] generates the 
deterministic part of the test particle dynamics, and \I/ determines the dissipative influence of 
the gas as an idealized noise reservoir for the evolution of the test particle. The Hamiltonian 
H is equal to jjP 2 + Hf, where Hf is an energy shift related to forward scattering with the 
gas reservoir particles and is defined below. 

The forms of \I/ and H, which I now define, are those suggested in [10] for a quantum linear 
Boltzmann equation modeling a test particle in dilute gas of distinguishable particles. The map 
\1/ is taken to be of the form 

= / dqe^T^e-^, (2.2) 

where are completely positive maps that act as multiplication in the momentum represen- 
tation 

T^p){p u p 2 ) = T ? (p' 1 ,p 2 )p(p 1 ,p 2 ), p e £i(L 2 (M 3 )), (2.3) 

for some functions T ? - : R 3 x I 3 C determined below. To define the maps T^, I must 
introduce some notation. Let m, M, = ™+ M refer to the mass of a single reservoir particle, 
the mass of the test particle, and the relative mass respectively. Let r\ be the density of the 
gas and r(k) = ( 2 ~^) 5e_5 '™ fc be the momentum distribution for a single reservoir particle. 
Define f : M + x [0 , it] — > C to be the scattering amplitude determined by the interaction 
potential between the test particle and a gas particle, and denote f(p2,Pi) '■= f(p, where 
P = HjPilh — IIP2II2 an d 9 is the angle between pi and p 2 . In my case, the scattering amplitude 
f(p,9) corresponds to the hard-sphere interaction, and I discuss its form in Sect. 12.21 Finally, 
the completely positive maps have a Kraus form 

T f (p) = f dvL^pLl,, (2.4) 

where (q)± is the orthogonal complement of the one-dimensional space spanned by q and the 
operators = L^(P) act as multiplication in the momentum representation as 

ill 1 

x ^x f + 2"'?). (2-5) 

V m M mm J 

In the above, the vector of operators P\\$ and P±$ have the forms P\\g = q and Pj_<? = P — P\\$. 
The Hamiltonian Hf will be given by 

H f = H f (P) = '- / dqr {q)Re [f{—q- -fP, — q- P . 2.6 

m* J R 3 m M m M ' J 
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Theorem 12. II is the main result of this article and states that the particle behaves diffusively 
for some diffusion constant D. A closed expression for D is discussed in Sect. 12.21 The conditions 
for the theorem are in terms of the initial density matrix p G B\ (L 2 (IR 3 )) . The statement 
that Y* pY 2 is trace class for unbounded operators Y 1 , Y 2 with dense domains D(YjJ, D(F 2 ) 
respectively, means that the corresponding quadratic form Q^i,^) = (^l^il/ 9 ^^) for ipi G 
D(Yi) and i\) 2 G 0(1^) is bounded and determines a trace class operator for its kernel. 

Theorem 2.1. Assume that Xjp, j = 1, 2, 3 are trace class and let pt(x) := p t (x,x) be the 
position density at time t. In the limit t — > 00, the rescaled density t^ p t {t^r), r G M 3 converges 
in distribution to a Gaussian with covariance matrix D I3, D > 0. Moreover, if in addition 
GiG 2 p and G\pG 2 are trace class for all Gi,G 2 G {X,-, Pj, j — 1, 2, 3} and |P| 3_n p|P| Tl is trace 
class for n — 0, 1, 2, 3, t/ien t/ie second moments satisfy 

t^ 1 / dx p t (x) XiXj = d~ijD + as t — > 00. 

2.2 Discussion 

2.2.1 Preliminary mathematical considerations 

The Lindblad equation (12. lft determines a trace-preserving collection of maps $t : i3i(L 2 (M 3 )) 
for t G M+ determined by $t(p) = Pi- The maps $t are completely positive, satisfy the 
semigroup property $£$ s = $t+ s , and are strongly continuous. The equation ( 12. ip has the 
generic form for the time derivative of a norm-continuous dynamical semigroup [31] except 
that the generator £ I consider is unbounded. It is therefore necessary to be cautious with 
respect to the mathematical interpretation of the Lindblad equation (12. ip . Analogously to 
unbounded master equations in classical probability theory [T3], a solution to an unbounded 
Lindblad equation may be constructible without being unique or preserving the trace the density 
matrices [6j. Intuitively, this would correspond to situations in which infinitely many jumps 
may occur in a finite time period and there may be "escape to infinity" in finite time. In the 
case of translation-invariant Lindblad equations, this problem can be reduced to a classical 
one [20]. The uniqueness of the solution is implied by the probability preservation for an 
embedded classical Markovian semigroup that governs the momentum distribution. In my 
case, the master equation governing the momentum density i*t{p) = Pt{p,p) has the form 

-rMp) = £o{vt){p) = / dpo(j(p,po)v t (po) - J{po,P)MP))i ( 2 - 7 ) 
at j R3 

where the jump rates J^(p+q,p) are written in terms of the differential cross section j^(p,Po) '■= 
\f(p,Po)\ 2 in the form 

m * 1 91 J(q)± dll\m M m M J 

X ''(" + 2 ^" + M^ (2 - 8) 

For the differential cross section ^ corresponding to a hard-sphere interaction, the so- 
lutions v t to (12. 7p converge exponentially fast in L 1 (M 3 ) to the stationary state v^p) = 
(2TrMf3~ 1 )~^e~2M^ 2 . This puts the dynamics (12. 7p in the "opposite" category as those in 
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which explosion occurs, since a particle starting with distribution ^o(-) = 5{- —p) for arbitrarily 
high momentum \p\ 3> Ma/3" must have distribution v t concentrated near a region of radius 
M2/3~2 from the origin after a finite time period. This is not surprising, since it holds for the 
classical linear Boltzmann dynamics modeling a test particle in a hard-sphere gas for which the 
marginal momentum distribution satisfies ( 12. 1\ with ^(p2,Pi) = 4:~ 1 a 2 for spheres with radius 
a. I return to a discussion of the existence and uniqueness of the semigroup $ 4 in the next 
section. 



2.2.2 Hard-sphere scattering 

Beyond the values m, M, rj, 0, the critical input determining the noise term \1/ is the scat- 
tering amplitude f(pi,P2) = f(p, &)• The scattering amplitude carries the essential two-body 
scattering information between the test particle and a single reservoir particle that becomes 
dominant in the low-density regime from which the linear Boltzmann dynamics arises (see [35J 
for a mathematical discussion of the scattering amplitude). The Hamiltonian for the two-body 
interaction is 

H 2 h 2 

Htot — ~ „ „ 7 A test — 7: A rcs + V{Xtest ~ ^res) 

2M 2m 

h 2 h 2 

-A cm -- — A re i + V(x Tel ), (2.9) 



( 00 


\x\ 


lo 


\x\ 



2{m + M) 2m, 

where V(x tes t — x res ) is the interaction potential, and A cm and A re i are the Laplacians with 
respect to the center-of-mass coordinates x cm = ^ x tes t + jj%res and x re \ = x res — x tes t- 
A hard-sphere potential of radius a is formally given by a potential 



V{x) 



Mathematically, a hard-sphere interaction is not defined as a perturbation of the V = case 
of (12. 9p . but as a certain self-adjoint extension of the Laplacian corresponding to the kinetic 
energy. For the region R = {\xt ex t — x res \ > a} = {\x re i\ > a} C M 6 , the Hamiltonian for the 
test particle and a single reservoir particle is 

h 2 h 2 h 2 h 2 

H tot = -^A tcst - — A rcs = - 2(m + M) A cm - ^-A rel (2.10) 

with Dirichlet boundary conditions on dR. 

Since the Hamiltonian (I2.10p is not a perturbation of a "free" Hamiltonian as in (12. 9p . 
the standard scattering formalism in [35] does not apply directly for the hard-sphere case. 
However, the scattering amplitude for a hard-sphere can be defined using the limits of the 
scattering amplitude for a permeable spherical potential V(x) = Vx(\x\ < a), V > 0. For the 
permeable sphere, the scattering amplitude has a partial-wave expansion 

ft 00 

f (P'^) = ^-E^ + ^^P) " 1 )^(cos(^)), (2.11) 
P £=0 
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where the values G C are determined by the equation below and £^ is the £th Legendre 

polynomial. The value Si(k) is a phase factor with 

^ 2) (k) log/fefV) - n\og'j e (xiK) 

0£{K) = rpj rpj , K G K + , 

^ log'/i^ >(k) - nlog' ji(nn) 

where n = (l -2m„Fa 2 /t -2 ^ 2 ) 2 , log' refers to the logarithmic derivative, and jg, hfp, and h!f^ 
are the spherical Bessel functions and Hankel functions of the first and second kind respectively. 
In the hard-sphere limit V — » oo, the above becomes 

St(K) = K G R+. 

h\ >{k) 

Even for an interaction as simple as the hard-sphere, there is no simple formula for f(p, 9) 
except when |p 1 and the much more challenging regime |p ^> 1. When |p <C 1, then 

f(p,0)«-a. 

In the regime |p ^> 1, there is a rough approximation [33] given by 

f(p,») « -\a(e-n**<) +i i±^)j l( £p sin ( S ))), (2.12) 

where Ji is a cylindrical Bessel function of the first kind. The left term e -2l t psm (§) is dominant 
except in a small region in the forward scattering direction 9 « 0. The second term in (I2.12p 
becomes dominant for 9 on the order of (|p) _1 <C 1 and yields a diffraction peak of height 
proportional to a(| p) ^> a and a series of much smaller peaks forming the classical Airy pattern, 
which vanish for 9 ;§> (|p) -1 . The total cross section <7 tot (p) = J n |f (p, 9)\ 2 , which characterizes 
the effective target area with respect to scattering, has the limiting values cr to t(p) ~ 47ra 2 and 
Ctot(p) ~ 2na 2 for |p 1 and |p ^> 1, respectively. The fact that cr to t(p) limits to twice the 
classical cross-sectional area ira 2 in the large p limit is known as the extinction paradox, since 
classical behavior would be expected in that regime (i.e. yielding the cross-sectional area na 2 ). 

The study of f(p, 9) through the partial- wave expansion when |p ^> 1 is well-known to 
be difficult, since the number of non- negligible terms to be summed in (12. lip grows on the 
order of |p. There is a successful approach based on exchanging the sum (I2.1ip for a complex 
integration over a contour for an analytic extension of the index I: 

h f dxx /ifdp) 



f(p,0) = / — - * Kh 7 £ A i(cos(0)). (2.13) 

; 2*p J c cos(TrA) ^)(| p ) ^ >> ^ > 

The contour C is is w-shaped around the positive real axis, although it must avoid enclosing 
any poles of h^\^p) as a complex function of A. The rewriting of the sum for the complex 



integration is called Watson's transformation [41 J. The complex integration formulation allows 
the application of the method of stationary phase and the calculus of residues in the evaluation 
of (I2.13p . This method was introduced by Watson and further developed in [TU [15]. Math- 
ematical detail can be found in j32[ [33], and [H] is an introductory book on the subject of 
scattering of waves from spherical objects. In addition to the quantum hard-sphere, the same 
formalism (12. lip applies for acoustical scattering and electromagnetic scattering from spherical 
targets. For the purpose of the current article, I essentially only require a small corollary of the 
above results, which is that the differential cross section ^(p,$) = |f(p,$)| 2 does not vanish 
away from the region of the forward direction 9 ~ as p — > 00. 
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2.2.3 The diffusion constant 



The diffusion constant D in Thm. l2TT1 is a sum of two parts D^ Q , Dj ps > that can be interpreted 
as driven by an averaged kinetic motion of the test particle and by spatial jumps made by the 
test particle. Let pi G M 3 be the classical Markov process with jump rates given by ( 12. 81) and 
starting in the stationary distribution v^. The constants D^ n and Dj ps have the form 

POO 

Ad* = 3- 1 / dtm Voo [x>(pt)-t>(po)], 

Jo 

Dj ps = 3~ 1 h 2 / dpu^p) I dqy^{di, n d 2 ,n^q){p,p), 

where di tn ,d 2tn , n = 1,2,3 are the partial derivatives with respect to the first and second 
arguments of T^(pi,p 2 ) respectively, and D : IR 3 — > M 3 has the form 

0(p) =(VH)(p) +u(p), for (2.14) 

U(P) =17 / ^Prel |Prel| r(— p rd + / ^ rel + |^ rel |^ Im [(d z f ) (p rel , fl)f (p rel , 0) . 

M J R 3 v m* M J n |p; el + \p iel \8\ L -I 

The integration J n is over unit vectors 6* G M 3 , is the angle between p re i and 0, and <9 Z is the 
mixed derivative d z = cos(2 -1 0)d p — 2p _1 sin(2 _1 0)d0. 

It is worth formulating classical analogies of the dynamics to better understand the contri- 
butions to the diffusion constant. Consider a classical dynamics in which the spatial variable 
X t G R 3 for the particle is driven ballistically by the momentum process p t through the disper- 
sion relation ( VH ) fpj : 



X t = X + [ dr(VH)(p r ) 
Jo 



The Markov process p t is time-reversible, since the jump rates (12. 8p satisfy the detailed balance 
property that the kernel A(pi,p 2 ) := <7(p 2 , p\)e^^~^ is symmetric in pi and p 2 . This follows, 
since 

Znm \p 2 — P\\ 

x / d y^(IIhLZ-2- 1 (p 2 -p 1 ), ^ v + 2-\p 2 -p 1 ))e~^ + ^ +p ^ 2 (2.15) 
J(p2-pi)± dQ\m m J 

and ^{pi-,p 2 ) is a function of the norm \pi\ = \p 2 \ and the angle between the vectors p\ and p 2 . 
Central limit theorems for integral functionals of time-reversible Markov processes were treated 
in [28]. A more general discussion on the subject of integral functionals of Markov processes can 
be found in the recent book [22] . This theory would suffice to show the renormalized processes 
N~^X]y t converge in law as Af -> oo to a Brownian motion B t in M 3 with diffusion constant 
D'L for D' > with the Green-Kubo form 



D' = 3-' dtE Uoo [( Vtf) (p t ) ■ (X7H) (p )] . 
Jo 
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Our expression for the diffusion constant D differs from the expression D' by the substitution 
of (Vi/)(p) with o(p) and the additional contribution D ]ps . Technically, the quantum dynamics 
is more closely related to a classical Markov process that moves with velocity (VH)(p\ and 
makes joint spatial/momentum jumps with Poisson rates that depend on the current state 
of the momentum (i.e. translation invariant rates). Suppose jumps (x,q) G M 3 x M 3 occur 
with Poisson rate density Rp~(x, q) when the current momentum is p G R 3 , and the following 
integration formulas hold: 

J(p + q,p)= I dxRp(x,q) and u(p) = / dx dq Rp(x, q) x. 

In other words, the marginal jump rates in momentum agree with the process p t and the 
averaged velocity generated by position jumps is u(p). Suppose also the rates Rp(x, q) satisfy 
rotation invariance Rp(x, q) = Rsp(Sx, Sq) for S G SO3 and the detailed balance condition 

Voo{p)Rp{x, q) = Voa{p + q)Rf + $(-x, -q). (2.16) 

Then, under some assumptions of the existence of moments for the spatial jumps, the diffusion 
matrix is D"I% for 



D" = 3- 1 / dtE Uoo [o(pt) • O(po)] + 3- 1 / dpv^p) / dxdqR p (x, q) 

Jo JR 3 JR 3 xR 3 



X 
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The detailed balance condition (I2.16P eliminates an extra term from appearing in the expression 
for D", and there is an analogous condition satisfied by the quantum dynamics considered in 
this article, which is discussed in Sect. 14.21 

Some similarity between Dj ps and the second term in the expression for D" can be seen by 
a comparison of Rp(x, q) with a Wigner transform of MfOpi , P2) ■ 

R'Jx,q) := (2tt)- 3 / dke^MJp- 2^hk,p + 2~ 1 hh. 

Although the classical analogy is not perfect, it gives some intuition for both the symmetries 
and complications in the analysis of the model. 

In the limit of large mass M for the test particle, the component of the diffusion 
constant dominates. Let A = |j, then as M — > 00 

^(^)* l^^^f ^ + *«T^) + 0(A3) ' 

where cr tot (p) is the total cross section defined above and cr z (p) = f n \d z f(p, 6)\ 2 . In the regime 
A< 1, -Djp S is small compared to the classical part D^, which has order A. 

I remark that the leading term for Dj ps in the limit of small A is not the same as the quantum 
diffusion coefficient D qdc discussed in [101 Sec. 5 ], which has the form 

c * = A f?(^)7 R / q<lV,ff °((T ) ' q) ' 
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where <7 (p) = J n (l — cos(0)) ff^(p, #)• The constant D qdc arises in a Brownian limit of the 
quantum linear Boltzmann dynamics ( 12. ip . The difference between the leading order of Dj ps 
and -D q( jc is a non- commutation between the diffusive and the Brownian limits. Most notably, 
the leading order for Dj ps depends on the derivatives of the scattering amplitude, but D qdc 
depends only on the values of the differential cross section ^(p ; ^) = | f (p>, 6>) | 2 . 

2.2.4 Remarks on the proof strategy 

The basic strategy I use for the proof of Thm. 12.11 is a standard one, which is formulated 
using a decomposition arising from the translation invariance of the model. The articles [261 
[3l [7] are similarly based on a translation symmetry. The translation invariance of the model 
implies a "decomposition" of the state space into fibers 23i(L 2 (M 3 )) « /® L^fl 

that evolve 

independently according to the dynamics. In other words, there are a family of maps p — > [p\t 
from i3i(L 2 (M 3 )) to L^M 3 ) for feel 3 and a family of semigroups $f : L^IR 3 ) such that 

The longterm diffusive behavior for the test particle is determined by information contained 
in the semigroups $^ for an infinitesimal neighborhood around k = 0. In particular, if 

~ (k) 

is the infinitesimal generator for & t , then the diffusion matrix is determined by the second 
derivative at zero for the eigenvalue e(k) G R_ + iR of having leading real part in the 
spectrum: D = — V® 2 e^|^ =Q . 

The analysis in the proof uses that the first few derivatives of in a neighborhood of k = 
are relatively bounded to Cq. The contribution to coming from the kinetic Hamiltonian term 
2^jP 2 is a multiplication operator with function gt(p) ■= jjk ■ p. The choice of hard-sphere 
interaction for the model has the advantage that the Hamiltonian term is relatively bounded 
to £q- Softer short-range interaction potentials will not have this property, and the argument 
would in the least have to be modified to cover other cases. 

3 Existence and conservativity of the dynamics 

In this section, I provide some technical discussion regarding the construction of the semigroup 
$t. The reader who is unconcerned by this issue may pass to the next section without any no- 
tations missed. A theory specific to translationally invariant Lindblad dynamics was developed 
in [20]. In Appendix |Aj I present a Fourier transform relating the noise map ^ to the maps T^- 
in ( 12. 2p . See [51 [21] and the review [22] for technical work regarding Lindblad equations with 
unbounded generators. 

I define the semigroup $ t : B\ (L 2 (1R 3 )) as the pre-adjoint maps corresponding to a semigroup 
of completely positive maps acting on the collection of bounded operators i3(L 2 (M 3 )) over 
the Hilbert space L 2 (IR 3 ) and satisfying the integral equation below: 

(m(G)^) = (<f>\G^)+ fdrC^, ^GD(P 2 ), G 6 £>(L 2 (M 3 )), (3.1) 

Jo 

^This equivalence is not meant to be precise. 
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where the form generator C : D(P 2 ) x B(L 2 (R 3 )) x D(P 2 ) — >• C is defined as 



C{4>- G; ij)) = (A<p\Gip) + (<f>\GAiff) + [ dq [ dv(e i ^L^(P\Ge i ^L^tP), 

JR 3 J(q) ± 

for A = -\H - 2" 1 ^*(I). The domain D(F 2 ) C L 2 (R 3 ) is the collection of all vectors 
whose density in the momentum representation satisfies J R3 dp\p\ 4 \(f){p)\ 2 < oo. The translation 
invariance of the dynamics is characterized in the form generator through the covariance relation 

C{(j)] e'^Ge^; ifi) = C(e Cbp (j); G; e* p ip), bet 3 . (3.2) 

This relation follows since A and are functions of P and by Weyl intertwining relations 
between e lbP and £k x . 

A semigroup $ t : £>i(L 2 (R 3 )) is said to be conservative if it preserves trace. This is equiv- 
alent to the adjoint maps sending the identity operator to itself (i.e. $^(1) = I), since 

Tr[$ f (p)] = Tr[^(I)p] = Tr[p], p G B 1 (L 2 (R 3 )). 

Conservativity is analogous to being stochastic for a Markovian semigroup [TBI Ch.X]. The 
weak*-continuity for $^ stated in Thm. 13.11 is equivalent to strong continuity for $ t , by a 
general result for semigroups on Banach spaces [3H Cor. 3. 18]. The translation covariance 
property (13.31) implies the same covariance for <3> t . 

Denote T m = {y E L 1 (M 3 ) | L 3 dp\p\ m \y(p)\ < oo}. In the proof of Thm. 13.11 I use results 
from Sect. [5] that the operator £ (defined in (12. 7j) ) is closed with domain 71 and generates an 
ergodic Markovian semigroup. The proof of Thm. [3TT1 reiterates some of the reasoning from [20J 
and concludes with an argument specific to my model verifying that Cq is equal to the generator 
of a certain Markovian semigroup defined through the quantum dynamics. 

Theorem 3.1. There exists a unique semigroup of completely positive maps $j : £>(L 2 (IR 3 )) 
satisfying the integral equation Iji3.1\) . The semigroup is weak*- continuous and conservative. 
Finally, the maps <3>j have the translation covariance property 

^(e-'^Ge^) =e- iip <5>* t (G)e iSP , beM 3 , (3.3) 

for allGe B(L 2 {R 3 )). 

Proof. By [6], a weak*-continuous semigroup of completely positive maps r £ < ^ nm ^* satisfy- 
ing (13. ip can be constructed if A generates a strongly continuous contraction semigroup on 
L 2 (M 3 ), and the domain D(v4) of A is contained in the domain of for all index values q, v. 
In my case, A is a perturbation of —jj^jP 2 by — iiJ/ — |^*(I). The operators ^*(I) = £(P) and 
Hf are both functions of the vector of momentum operators and are relatively bounded to \P\ 
(see Lemmas l5.1l and lC.ll) . Thus, A is m-accretive with domain D(P 2 ) and generates a contrac- 
tive semigroup. Finally, the operators are relatively bounded to \P\ by Lem. IC.ll and thus 
include D(P 2 ) in their domains. It follows by |6J, there exists a solution <jj mm )'* to (13. ip such 
that for any second weak*-continuous collection of completely positive maps %(G) solving (13. ip 
that 

%l min) '*(G) <% t (G) <\\G\\I, GeB(L 2 (R 3 )), G > 0. 
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The semigroup <xj mm ' ) '* is referred to as the minimal solution. From [20], the relation ( 13. 2 j) 
has the consequence that the minimal solution <j( mm )'* satisfies the translation covariance rela- 
tion ( 13~3|) . 

With the construction of ^( mm )'* ; the next question is of its uniqueness as a solution to (13. ip . 
By [6], the uniqueness of the solution <j( mm )'* j s implied by *(I) = I f or all t > (i.e. 

^(min),* . g conservative). As observed in [20], the translational covariance of <£( mm )>* reduces 
the problem of checking conservativity for the Lindblad dynamics to checking stochasticity 
for a classical Kolmogorov semigroup. I prove below that the Kolmogorov semigroup is e tc o. 
By Prop. 15.31 e * £ ° is ergodic to the stationary state (27vM(3~ 1 )~2e~2M^ 2 , and it follows that 
e tc o(l R s) = 1 R 3. This would entail that := <j( mm )'* j s the unique weak*-continuous solution 
to (13~TT) . 

The algebra A = {g(P) \ g G L°°(R 3 )} is invariant under the action of < £( mm )'* ) since bounded 
functions of P are the only elements in i3(L 2 (M 3 )) that commute with the Weyl operators e lbP 
for all b G M 3 . Hence, I can define a semigroup of maps : L°°(R 3 ) through g t = S£(g) for 

g t (P):=^r n ^{9(P)), geL°°(R 3 ). 

The semigroup <S£ is contractive and maps positive functions to positive functions. Moreover, 
^(Ir 3 ) = 1r 3 if and only if x£ ml n ^'*(I) = I. Let L* be the generator for S% and L be the adjoint 
generating S t . Denote the dense domain of L by D(L) C L 1 (R 3 ). I will show that L = £ - 
By differentiating flU]) at t = for G = g(P) with # G L°°(R 3 ), 



It 



t=o 



dpy(p)S* t (g)(p) 

dq I dpy(p)J{p + q,p) g(p + q) - dpy(p) dqj(p + q,p)g(p) 
dpC (y)(p)g(p), (3.4) 



where y(p) = <fi(p)ij)(p). The set of possible y(p) with </>, ip G D(P 2 ) is 71. As a consequence 
of Part (2) of Lem. 15.31 the generator Co and the multiplication operator (Wly)(p) = \p\y(p) 
(having domain 7i) are relatively bounded to each other. In other words, there is a c > such 
that for all y G 7i the graph norms for £ an d 971 satisfy 

c-^IMIi + < \\y\U + ||£ (2/)||i < cQli/ll! + ||SDt(y)||i). (3.5) 

This implies the right side of (13 ,4p is bounded for y G 7i- Moreover, any closed operator 
agreeing with £ on a subset of 7i that is dense in 7i with respect to the graph norm of OK 
must be equal to jC q . Continuing with (13 ,4p . for all y G 71, g G L°°(IR 3 ) 



=0 



dpy(p)S*(g)(p) = [ dpS t (y)(p)g(p) 

dt t=o J R 3 

dpL(y)(p)g(p) (3.6) 



Hence, L*(y) is the weak limit of t 1 (S^(y) — y) as t — » for all y E T^. Consequently, by the 
remark [27, IV. 1.5], 71 must be contained in the domain D(L). Moreover (13 ,4p and (I3.6P imply 
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L agrees with £ for elements in %. Since L generates a contractive semigroup, it must be 
closed. Finally, 74 is dense in 7i with respect to the graph norm of 971, and therefore £q = L 
by the remark above. 

□ 

4 The fiber decomposition and detailed balance 

In this section, I discuss a few symmetries and decompositions for the dynamics. The dynamics 
also satisfies a rotation covariance, which I do not explain. 

4.1 Fiber maps 

The Lindblad generator C satisfies translation covariance with respect to position shifts a G M 3 

C(e-^ p pe^ p ) = e< p C{p)e^ p . (4.1) 

A consequence of this symmetry is that the dynamics admits a fiber decomposition in which 
individual fibers obey autonomous dynamics. Similar fiber decompositions can, for instance, be 
found in the study of periodic Schrodinger potentials [361 Ch. XIII. 16] and in the study of certain 
translation invariant noisy dynamics for a quantum particle living on Z d , [261 El U\- Also, the 
autonomous Markovian evolution for the momentum distribution (12. 7p is a special case of the 
fiber decomposition and provides a useful tool for proving a solution of an unbounded Lindblad 
equation is unique when the Lindbladian is translation-covariant [20] (as discussed in the previ- 
ous section). In this and future sections, I will denote T n = {y G L X (1R 3 ) | J R3 dp\p\ n \y(p)\ < oo} 
and use the following notations for norms 

• || y ||p is the standard L P (1R 3 ) norm for p G [1, oo]. 

• ||p||i := Tr[|p|] is the trace norm for elements in B\ (L 2 (IR 3 )) . 

• || A || is the operator norm for elements A G B(L P (R 3 )). 

• \\y\\A is the graph norm for a densely defined operator A on L 1 (M 3 ): 

\\y\\ A := + 



Now, I define the fiber maps. For k G M. 3 , there is a contractive map p — > [p\t that sends 
an element p G i3i(L 2 (IR 3 )) to an element L X (IR 3 ) satisfying the formal relation 

[PUP) = Pip- P + 2^)' ( 4 - 2 ) 

where the kernel of p on the right side is in the momentum representation. A rigorous definition 
of [p\r requires something more than the integral kernel of p, which is only defined a.e. 1R 3 x 
M 3 . However, I can define the contractions : £>i(L 2 (M 3 )) — > L 1 (1R 3 ) through the Riesz 
representation theorem. Let g G L°°(IR 3 ), then by standard properties of the trace 

\Tr[^g(P + 2- 1 hk)p} \ < \\g(P)\\ \\p\U 

= IMUMIi. (4.3) 
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By the Riesz representation theorem, there is a unique complex- valued measure p^ on R 3 with 
total variation < \\p\\i and satisfying the first equality below for all g G L°°(R 3 ). 

Tr[e rk}l g(P + 2- 1 hk)p}= [ p %p (dp)g(p) 

dp[p}k(P) 9(P) (4-4) 

The continuous spectrum of the family P^, j = 1,2,3 implies the measure p^ must be abso- 
lutely continuous with respect to Lebesgue measure, and I denote its Radon- Nikodym derivative 
by [p\t G L 1 (R 3 ) in (14.41) . It is clear from (14. 2 p the functions [p]^, k G R 3 also determine the 
density matrix p. 

The fiber decomposition for the dynamics $ t can be characterized by the existence of semi- 
groups $p : L : (R 3 ) such that 

[*«(P)]* = *?W (4-5) 

The relation H4.5f) is equivalent to noticing that maps the Banach space Bj: = {e lkX g(P), g G 
L°°(R 3 )} C i3oo(L 2 (R 3 )) into itself for each k with 

^* t (e ik)l g(P + 2- 1 hk)) =§f ) '\g)(P + 2- 1 hk). (4.6) 

where $[ fc - ) '* : L°°(R 3 ) is the adjoint semigroup of $^ . 

Lemma 4.1. For a// G R 3 ; t/ie semigroups $^ : L 1 (R 3 ) are strongly continuous and con- 

— (k) 

tractive. Consequently, the generator £g of the semigroup <5> t is a closed operator with dense 
domain D(£g) C L X (R 3 ) ; and all z G C with Ke(z) > belong to the resolvent set of with 

~ (k) * 

Proof. It is equivalent to establish the contraction property for the adjoint semigroup $J . 
The supremum norm of an element / G L°°(R) is equal to the operator norm of f(P) G 
S(L 2 (R 3 )), so I have the first equality below for g G L°°(R): 

ll$f ) '*«l|oc=||$P(0)(P)|| 

= \\^{e' M g(P + 2- l hk))\\ 
< ||o( J P + 2- 1 M)|| = \\g\U 

The second equality above is by the defining relation ( 14. 6p . and the inequality is due to 
being a contraction. By similar reasoning, the weak*-continuity of §* t implies weak*-continuity 

for , and consequently the pre-adjoint semigroup $^ is strongly continuous. 

The consequences for the generator and its resolvents are standard and can be found 
in [23 Sec.IX.2-3] or [TJl Sec.XIII.9-10]. 

□ 
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By the same reasoning as for Co in the proof of Thm. 13.11 the domain D(£g) must include 
Ti- An implication of Prop. 15.31 is that the domain contains T\ for all k. The operator Cj: is a 
sum of terms 

C n =~h n + J n -S n (4.8) 

where the terms J7g, £t : 7i — > L 2 (R. 3 ) arise from the Hamiltonian, and gain and loss parts 
of the Lindblad generator ( 12. ip . respectively. The terms above are determined analogously to 

$f } in ((431) such that for p in a dense domain of i3i(L 2 (lR 3 )), 

Hp\% = l[BMh> m% = [*(p)] 5 . E M% = tr^mp}]*- 

Formally, I have 

h % (f)(p) = (H(p-2- 1 hk) - H(p + 2- 1 Hk))f(p) 

W){P) = I dp Jz(p,p )f(p ) (4.9) 

£ % {f){p) = 2- 1 (£(p-2- 1 hk)+£(p + 2~ 1 hk))f(p) 

where £{p) = L 3 dpo JoipoiP), and Jj:(p,Po) is complex- valued pseudo jump kernel having the 
form 

Jk&Po) = M^ Po (p - 2- l hk, p + 2~ 1 hk). 
The jump "rates" J%{p,Pq) can be written more explicitly as 

= Mf(p- 2- l hk, p + 2~ l hk) 

= / dvL^(p-2- 1 hk)L^(p + 2- 1 hk) 

J (a) , 



where Q± and 6± are defined as 



^kT 1 e-^t / d«rr(tr+ 2- 1 -? + ^P| k -)f(Q_, e_) f (q + , e + ) (4.10) 



Q ± = Q ± (g,^p,fc)=((^-^fc±2- 1 ^)) 2 + 4-V) : 

e 



U M 



2 — 1 H 

± = Q±(q,v,P,k) = 2 tan" 1 f- ! ^ ^ -) 



«± = z±(q,v,p,k) = v + 2 _1 — <f+ — (%±2~ /*%), 

771/^ lvJ. 



and I have used the simple relation 



\I l-> \- 13 „ k 2 /-> -1 m -> m - \ 

r(z_)2r(z + ) 2 =e s« 2 u*rlv + 2 — q H pii^). (4.11 1 
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4.2 Detailed balance 

There is also a detailed balance symmetry for the noise map The function T^(pi,p 2 ) for the 
multiplication maps : ,Bi(L 2 (R 3 )) (defined in d2.4[) ) can be written in the form 

1 _L 

T f (pi,p 2 ) = 1 1 7<APi,P2), (4.12) 

where J^ipi,^) is the kernel for a positive operator on the Hilbert space L 2 (R 3 ) for each q and 
satisfies the symmetry 

7?(Pi,p 2 ) = J-q(Pi + q,P2 + q)- (4.13) 

The explicit form for 7g(pi,j?2) can easily be worked out, and along the diagonal, it reduces to 
IqiPiP) — A(p + q,p) for the kernel A defined in (12.1 5p . 

The above translates into a decomposition for the noise map Let T : i3i(L 2 (M 3 )) by 

given by T(p) := v£o(P) pv£o(P), and formally define T = T _1 \l/T. The above decomposition 
for the maps implies that T satisfies 

f^) = r(p), peB^iR 3 )), 

where the bars refer to complex conjugation is in the momentum representation. We do use 
this formal relation, although the symmetry (14. 13|) is useful for seeing a cancellation in the 
proof of Prop. 17.11 

5 The "marginal" momentum dynamics 

For a translation-invariant linear Boltzmann dynamics, the momentum is a Markovian process. 
Thus, the marginal probability density for the momentum obtained by integrating out the 
spacial degrees of freedom is governed by an autonomous master equation. The quantum 
situation is similar in this respect as seen in the classical Kolmogorov equation (12. 7p having 
jump kernel J^{p2,Pi), which corresponds to the fiber decomposition from the last section for 
the k = case (I drop the index for k — 0). I denote the gain and loss terms for the Markov 
generator Co by J and S, respectively. In this section, I will characterize the ergodicity for 
the classical stochastic semigroup e tc ° and collect some analytic facts about the generator Cq. 
As mentioned in Sect. I2.2[ the momentum distribution is exponentially ergodic to a Gaussian 
with variance and that statement is made more precise below in Prop. 15.31 In the analysis, 
I rely on soft consequences of Nussenzveig's results in [33] to ensure the scattering cross section 
^(p, 6) does not vanish for 6 e (0, 7r) bounded away from the endpoints as p /* oo. 

The following lemma gives upper and lower bounds for the escape rates £(p) = L 3 dp J'(p ,p). 
An implication is that the operator \I/*(I) = £(P) acting on the Hilbert space L 2 (M 3 ) is rela- 
tively bounded to \P\. As I mentioned in the introduction, the total cross section cr tot (p) > 
appearing in the bounds of Lem. 15. II approaches Aria 2 for p 1 and 27ro 2 for p ^> 1. 

Lemma 5.1. The escape rate £{p) has the following upper and lower bounds 

£_V(£\p\)<£(p)<e + ?\p\, 
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where the constants £, S.^J 7 are defined as 

8 = ( sup a tot (p)) — (— ) (-— ) 2 , J-= ( sup a tot (p) — — . 
P eM+ m * m pgR + M m 

Proof. The escape rate can be written 




where p re \ = — v — jjp±q + 2 1 q and 6 = 2 tan 1 ( , m „ ^ ifU — r) . The integral J Q is over all 

unit vectors S 6 K 3 and 6 is the angle between p re \ and The second and third equalities 
in (15. ip are changes of integration variables, where the third equality is simply due to p re \ — > Po = 
^-Pmi+jjP- For the change of integration in the second inequality, I consider a smooth extension 
of ^(|Prei|,0) to a function ^(p re i, Prei.o), Prei.o e M 3 , where ^(|p rel |,6») = f^(p re i, Prei.o) when 
IPreil = |Prei,o|- I relate the variables p re i i0 and p rei to v, q e M. 3 as p reliQ = ^fv - ^p±$- 2~ 1 g 
and as above for p re \. I can go from integration over p Teh 9 to integration over p re i,Prei,o an d 
then to integration over q, v by the same reasoning as in (101 App.A]: 




With ( 15. ip . I have the following upper bound for the quantum escape rates: 
£(p) < ( sup or tot (p))ri— / dp \ — po - -rjP\— — -y 

<( JS pa t- (p)),^(2-(^). + -H). 
This gives the coefficients for an upper bound. For a lower bound, I have 

£ (PJ > ( inf or tot (p))i/— / dpo — Po-^jP\— ~ 7T^"' 

The right side is minimized for p = with the value £. Moreover, since the integral is a 
weighted average with mean zero of the convex function F(po) = r^fjpb — 3fPl> the integral is 
smaller than the value -F(O) = jj-\p\- 

□ 
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The following theorem summarizes a small part of the results from [33]. To prove Co has 
a spectral gap in Prop. 15.31 it is imperative that the scattering does not all become absorbed 
in the forward direction when the test particle is traveling with high momentum. The lemma 
states that the differential cross section ^(p, 0) approaches 4 _1 a 2 as p — > oo for 6 bounded 
away from the endpoints 0, n. 

Theorem 5.2 (Nussenzveig). As |p — > oo, there is uniform convergence 

£ = |,(p,,)|»_4 ,<9<r-6, 

for any fixed 5 > 0. 

Proof. By [33l Sec. 3], there is uniform convergence 

f( P) 0) _^ _^ e -2if PS m(f) 

as |p — > oo for 9 bounded away from the endpoints of the interval [0, it]. The estimates in [33] 
even allow for uniform convergence over intervals [(f p) -7 , 7r — (f p)~ 7 ] for < 7 < 3 _1 . 

□ 

For y G L 1 (1R 3 ), define the multiplication operator $Jl(y)(p) = \p\y(p). Part (5) of the 
proposition below states that the graph norm for 9Jt n is equivalent to the graph norm for the 
operator Cq. This equivalence is useful, since || • ||ajp. is easy to compute in some situations 
and the norm || ■ H^n is contractive under the dynamical evolution e tc °. By the discussion in 
Sect. El the domain D(£ ) contains %, and the domain for D(£^) also includes % by the same 
argument. Since £g is closed, Part (5) of Prop. 15.31 for n = 1 implies D(£ ) — 7i, and Part (6) 
implies D(£^) contains 7i- Moreover, Part (5) also implies £ n maps elements in Tn into 7/v-n 
for N > n. 

Proposition 5.3. Let £g be the generator for the contractive semigroup $j acting on the 
k-fiber. Define the projection P to operate as Py = ( L 3 dpy(p))u 00 for y E L 1 (IR 3 ). 

1. Zero is a non- degenerate eigenvalue for jCq with corresponding spectral projection P. 

2. The remainder of spectrum satisfies XX^o) — {0} C (—00, — g] + iR for some q > 0. 

3. Given 5 > 0, there is a C 5 > such that ||$| 0) (I - P)|| < C a eW for allt>0. 
4- For C$ > defined above and z G C with Re(z) > 5 — g, 

11 1 11 ^ 1 1-1 , ^ j 11 ^° 11 ^ o 1 

< U + - ; ; and — < 2 1 



11 C — z" Re(z) — 5 + g 11 C — z" Re(z) — S + g 

5. The graph norms || ■ H^n and || • ||gjt« are equivalent. 

6. Cj: is relatively bounded to Cq- 
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Proof. 
Part (1): 



It can be verified that the state Voo{v) = (27rM/3 _1 ) 2 e"2M^ satisfies Cqu^ = by a direct 
computation. This also follows from the detailed balance form of the jump rates (12.151) . Since 
the jump rates J {02, Pi) are strictly positive, the process is ergodic and the only stationary 
state is z/qo. Thus, is a non-degenerate eigenvalue of Co. 

Part (2): 

Let I r : L X (]R 3 ) — > L X (]R 3 ) be the projection onto the ball of radius r > (i.e. multiplication 
by the indicator function x(\p\ — r )) an d define the operator Q r = £o(I — I r ) — £ I r with domain 
T\. I will argue the following points: 

(i) . J\ r is compact. 

(ii) . For sufficiently large r, Q r has spectrum in (— oo, —h] + iR for some h > 0. In fact, 

\\e- tGr \\ <2e- th , (5.3) 
where || • || is the operator norm for maps on L 1 (IR 3 ). 



Assuming the above two points, then since C = Q r + J7"l r , it follows by Weyl's theorem 
that C and J\ r have the same essential spectrum. Since £ generates a contractive semigroup 
XX^o) C — R + + iR, the set XX^o) H [(— h, 0] + iM] may only contain elements in the point 
spectrum of Co and no accumulation points. In fact, there will only be finitely many elements 
in XX^o) H [(— h + e, 0] + iM] for any e > 0. To see this, let us suppose there were an infinite 
sequence of normalized eigenvectors e n G L 1 (M 3 ) for C with eigenvalues v n G (— h + e, 0] + iM. 
Since there are no accumulation points, the v n must have |Im(i; n )| — > oo. Notice that 

Jl r e n = (C - Q r )e n = (v n - Q T )e n . (5.4) 

The left side of (15. 4p tends to zero in L X (IR 3 ) as n — > oo. To see this, first observe 

||i7Ire n || i< ( sup ll^lli) / \e n {p)\, 

V |pl<r J J\p\<r 

where lp (p) = J{p,Po)- By using that Coe n = v n e n and ||e n ||i = 1, there is the inequality 

sup \e n [p) \< — — — ^0 as Im(w„) — > oo. 

|p|<r 

Moreover, the right side of (15 .4p must have L 1 norm > 2~ x e for each n, since (y n — Q r ) 1 = 
J Q °° dte t ^ r ~ Vn ^ has operator norm < 2e _1 by (15. 3j) . This contradiction shows there can be 
only finite many eigenvalue in the region (—h + e, 0] + i R. All of these eigenvalues must have 
real part strictly smaller than by the mixing generated from the strictly positive jump rates 
J{j>2iP\) > 0. Hence, there is a gap q between and the real component of the remainder of 
the spectral values £(£o) ~ {0}. 
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I now turn to establishing the points (i) and (ii). For the compactness of Jl r , notice Jl r 
maps elements y G L^IR 3 ) to integral combinations of the functions lp~(po) through the formula 



Jl r (y) = / dpy(p) l p . 

J\v\<r 



'\p\<r 

The compactness of Jl r follows, since lp G L^R 3 ) is a continuous function with respect to the 
L 1 norm over the compact set \po\ < r. 

Proving Re^((? r ) is bounded away from zero for large enough r is more involved. Notice 
that the semigroup e t<3r : L 1 (M 3 ) maps positive functions to positive functions, since Q r agrees 
with £ on the image of I — I r and acts as a real- valued multiplication operator on the image of 
I r . Showing that y t = e tSr y has \\yt\\i < 2e" i/l ||?/||i for h > and all positive functions y gives 
the upper bound Re(^ Q T ) < —h. 

It is convenient to consider a slightly deformed norm || • || ^ that is equivalent to the ordinary 
|| • ||i norm. Define the function w(p) = 2 — -^-j- and the corresponding weighted norm 



\y\\ w = [ dpw(p)\y(p)\, yeL l (R 3 ), 
Jr 3 



so that I have ||y||i < \\y\\ w < 2||y||i. For a positive element y G 7i, 
d 



dr 



yt\\ w = / dpw{p) (G r yt)(p) 
Jm. 3 

dp(G*w)(p)y t (p) 
dp(£* w)(p)y t (p) - I dpw(p)£(p)y t ( 



_]>r J\p]<r 

Since S(p) is bounded from below by £(0) = £_, I only need to demonstrate sufficient decay 
from the first term on the right-hand side. My estimates will show there are h',r > such that 

(C* w)(p) = [ dp (w(po)J(po,p) - w(p)J(po,p)) < -ti V(|pl > r). (5.5) 
Jm. 3 

I can then take h — h! A £_ as a lower bound for the rate of exponential decay for \\yt\\ w an d 
thus also for ||^||i. The physical reason for expecting the exponential decay in \\yt\\ w is that 
when the test particle is traveling with high enough momentum p , then the first collision will 
typically knock it to a momentum pi such that is a fraction of the value \po\. How large po 
should be for this effect will depend on the masses m, M and the temperature (3, in order that 
the test particle is traveling at a much higher speed than the particles in the reservoir. Finally, 
this effect is amplified by the fact that collisions occur more frequently in proportion to the 
speed of the test particle due to the "hard" character of the interaction. The combination of 
these two features is enough to be visible to the gently sloping weight w(p). 



/ dp (w(po)J(p ,p) - w(p)J(p ,p)) 
Jr 3 

m f dq f ,-,/., n _i _ m _ x 

=r) ^ / i^i / dvr ( v + 2 —y + M p w 

m * Jr 3 \Q\ J(q) ± m * M 



X \f{\Prel(q,V,p±q)\,9(q,V,p ± j))\ 2 (w(p + $ - w(p)) 
: / / dp rd \p rel \r( p rel + — p) \f(\Prel\,0)\ 2 (w(P + Prel ~ 0\Prel\) ~ w(pj) , (5.6) 



19 



where 9 is a vector variable on the unit sphere with angle 9 from p Ie i, and the second equality 
above follows by the same change of integration as in the proof of Lem. 15.11 

By changing variables, b = ^p r ei + f| V an d rewriting the integrand from (15.61) . I have 



Due to the Gaussian factor r(b), the integrand decays rapidly for |6| 3> (j) 2 - Hence for large 
\p\, I have \p\ 3> \b\ and 



mv-^b\{\v\-\ 






- 2^61 1) 


(b1 + i)(| 


{l-^)p+^b-9\ 





I note that 



(1 — — H — — cos(6>)) + — £ sin 2 = (1 - — ) 2 + + 2— (1 - — ) cos(#) 

_ m 2 + M 2 + 2mM cos(0) 
(m + M) 2 

<1, 

where equality occurs only for 9 = 0. 
With the estimate above 

f -> TYV* [ -* -* 

I dp (w(po)J(po,p) - w{k)J(p ,p)) = -V^tj / dbr{b) 
Jm3 m 6 M J r3 

x / |f(|™ 6 --^, e )P( !H±>1 r - l) + OdpT 1 ). (5.7) 

■/n m * M v (m 2 + M 2 + 2mMcos(#)) 5 y 

— * — * 

Since r{b) is probability distribution concentrated around b = 0, I have an average of terms of 
the form 

/ |f(p,#)| 2 ( r-l) for | P |»1. (5.8) 

Jn v (m 2 + M 2 + 2mMcos(#)) 5 J 

However, by Thm. I5.2[ ^(p,^)! 2 converges uniformly to 4~ 1 a 2 as p — > oo for 9 bounded away 
from 0,7r. Hence (15. 8 j) will include a contribution from the integration of 9 away from the 
boundary points of approximately 



a 



2 f ( m + M \ w .mAM 



M ( (m+M) ^r- 1 ) 



4 7o V ^ m 2 + M 2 + 2m M cos(0)) 5 ' V mM 

,ra A M 



2 ' 

7ra 



m V M 

Hence, there exists an r > and an h! > such that (I5.5P holds, which completes the proof. 
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Part (3): 

Since the projection P commutes with the maps the operators — P) form a 

semigroup that may be viewed as operating on the Banach space Im(I — P) C L 1 (M 3 ). By J27J 
Sec. IX. 4], it is enough to show there is a C$ > such that 

iK^n'-^N^ih (5 - 9) 

for all z > —q + 5. By Lem. I4.1[ for z > I have the second inequality below 

\\(C - zy\l-P)\\ < ||(£ -^) _1 || +^ 1 ||P|| < 2z'\ (5.10) 

On the other hand, the operators (Co — z) 1 (l — P) have finite norm for z G (— fl, oo), and the 
norm depends continuously on the parameter z. The continuity gives 



max 

-g+S<z<l 



(C -z) \l-P)\\ <oc, 
and by combining this with (I5.10p . I can find a C$ so f)5.9p holds. 
Part (4): 

This follows by Part (3), the fact ||P|| = 1, and the formula 

POO 

(Co -z)-\l-P) = - dte- tz <bf\\ - P) 

Jo 

for Re(z) > 5 — g. 
Part (5): 

To bound || • \\c using || - Hon, I observe 

\\£ov\\i < 2||fy||i < 2(F||2/||i + ^|| \p\y\\i) < 2(S V 7) \\y\\ m , (5.11) 

where the first inequality uses ||i7y||i < ||£y||i (with equality when y > 0), and the second uses 
Lem. 15.11 Hence, ||£ y||i is smaller than a constant multiple c = 2(8 V J 7 ) of Hi/Hot- To extend 
to n > 1, I will use that there are C/s such that for all y G L 1 (M 3 ) 

\\[^,Co]y\\i<C j \\y\\ mj+ i. (5.12) 

This can be seen by the inequalities 

\\[m j ,Co}y\U 

< 



dpy(p) / dqj(p + q,p)(\p + q\ J - \p\ 3 ) 

/ / rf^rel |Prel| / (|p + Prel - 0|Prcl| \ 3 - \p\ J )r( p rei + —p) — ( \p xe i \ , 9) 

J R3 J R 3 J n v m» M J dil 

<2 3 — sup otot(p) / dp\y(p)\ / dp Te[ \p reX \{\p\ 3 + 2 J \p rel \ J )r( — Prei + jjp), 
m* P eM+ Jrs J R 3 m* M 
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where the equality uses the integral definition (12.81) of the jump rates J(p2,Pi) and a change of 
integration from q, v to variables 9 € Q, p rc \ 6 i 3 , as in the proof of Lem. 15.11 where 9 is angle 
between 9 and p re i. The second inequality simply uses (a + 6)- 7 < 2^(a^ + hP). By performing the 
Gaussian integration, the above is bounded by a linear combination of 1 1 £DT*2/ 1 1 1 for i < j + 1. 
By interpolation || £o]?/||i is bounded by a constant multiple Cj of ||y||a>tM-i. 
Using (15. lip and ( 15. 12ft . ||£o2/||i is bounded by 

Ky\\i < c\\C™- l y\\ m 

< c\\Ct l y\\i + c\\Co WlC^yh + c\\[C , ^]Q~ 2 yh 

< c\\Q- 2 y\\ m + (c 2 + cCJWC^vWw- (5-13) 

With interpolation to bound || ■ by || • ||ot 2 ; I can then proceed by induction until the right 
side is 

Now I must to show || ■ ||gjt™ is bounded by a constant multiple of || ■ \\cn. The result will 
follow by an induction argument if I can show 

II \p\ n y\\i < CJ bT^lli 

< C' n (\\ l^^lllblli+ll Ip-T^lli), (5-14) 

where Q T and J\ r are defined as in Part (2), and the first inequality holds for some C' n > 
for large enough r. The second inequality is the triangle inequality using Q r = C — J\ r) and 
the operator norm || Ipl™ -1 ^!,.!! is finite by the calculation below. Since \p\ n Jl r has a positive 
integral kernel, its operator norm as a map on L 1 (M 3 ) is bounded by the following 



|| \p\ n Jl r \\ = sup / dp- \ph\ n J(p ,p) 

\p\<r JR 3 

= SUp / dp le i / \p + p rc i - 9 % cY | I "r( p re[ + — p) — ( \p Te y \ , 9) 

m* \p\<rJiR3 Jq v m* M dll 

71 I Til TYl 

<2 n — sup a tot (p) sup / dp rel (\p\ n + 2 n |p rel | n )r( — p rel + —p) < oo, 
peR+ |pl<rJR3 v m* M 

where the second equality is the standard change of variables. The inequality is from |a + 6| n < 
2 n (|a| n + |6| n ) for a = p and b = 9\p TC \\ —p TC \ along with the fact \9\ = 1. The right side is finite 
since |p| < r and the integrand is a Gaussian against a polynomial. 

To show (I5.14p . I begin with an application of Lem. 15. II to get the first inequality below 

HbTylh < ZT'W W'Syh < ir'W IpT^IU 

< (KD^W W~ x Grv\U < 2(h n Z)- 1 \\ \p\ n - l g r y\\i. (5.15) 

The second and fourth inequalities follow from || • ||i < || < 2|| • ||i. The critical input in ( 15. 1 5f) 
is the third inequality, which follows from the triangle inequality and similar reasoning as for 
the proof of (15. 5p to get the first and second inequalities below respectively, 

ii w^Qrvh > ii ipt^iu - ii IpT~V(i - i^iu 

> h n \\ \p\ n - l £y\\ w , (5-16) 

for some h n > 0. I then set C' n := 2(h n! F_)~ 1 . The use of the norm || • \\ w is only required for the 
n = 1 case to get the second inequality of (I5.16p . and the physical motivation for the second 
inequality is the same as for (15.51) . 
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Using (15.14p inductively, then ||y||aji" w iU be bounded by a multiple of ||?/||,c«. 
Part (6): 

To show || • \\c-. is smaller than a constant multiple of || • ||£ , by Part (3) it is sufficient to 
prove ||£j$/||i < c lblbt for some c > 0. By (14. 8ft . Cj: is sum of a Hamiltonian term —lh a 
gain term J7g, and a loss term £t- The Hamiltonian term is multiplication by 

hiP) = ~k ■ p+ {H f (p- 2~ 1 hk) - H f (p + 2- 1 hk)) . 

The first term is linear in p and Hf(p) has a linear bound in \p\ by Lem. IC.ll The loss term 8$ 
is multiplication by 

2- l £(p-2- 1 hk) + 2- l £(p + 2- l hk), 

which is linearly bounded in \p\ by Lem. 15 . 11 For the gain term the values J%{p + q, p) are 
bounded by 

| MP + Q,P)\< ( dv\ L^(p- 2~ x fJz) Lq,v(p + 2~ 1 /^) \ 

< 2- 1 [ dv(\L^{p- 2~ 1 hk)\ 2 + \L^{p + 2~ 1 hk)\ 2 ) 

= 2- 1 (J(p- 2- x fik + q,p- 2- l hk) + J(p + 2~ 1 hk + q,p + 2~ 1 hk)) . 

The above implies 

WJkVh = / d P dq\j(p + q,p)\\y(p)\ 

< 2' 1 I dp(£{p-2- 1 hk) + £{p + 2~ 1 hk))\y(p)\ 

< I dp(~8 + 2~ 1 T\p-2~ 1 hk\+2~ 1 T\p + 2- 1 hk\)\y(p)\ 

< {£ + 2- x T\k\)\\y\\ x +T\\m y \\ x , 

and the right side is smaller than a multiple of ||j/||gn- 

□ 



6 Smoothness of the fiber generators 

The main goal of this section is to show the operators £g have three derivatives in k that 
are relatively bounded to Cq. Showing this will require bounding derivatives of the scattering 
amplitude (<9^f)(p, 9) for N < 3 for large and small p. By the discussion in Sect. 12.21 f(p, 9) 
has a peak with height proportional to p occurring in the forward direction 9 <C 1 for p ^> -. 
Hence, it would not be expected that the supremum sup e6 j 0j7r ] |(c^f)(p, 9)\ is well-behaved for 
large p. However, controllable behavior can be found by integrating the square modulus over 
the angular variables (dQ = sin(9)d9d<f)) 

[ |(aff)(p ; 0)| 2 = 27r [* d9sm(9)\(d?f)(p,9)\ 2 . 
Jn Jo 
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The analysis in the proof of the lemma below boils down to standard relations for Bessel 
functions and Legendre polynomials. 

Lemma 6.1. There is a unitless constant C > such that for N = 1, 2, 3 ; 
1. 



p> 

2. 



sup / \(d^)( P ,9)\ 2 <a^rC, 



snp / |(^f)(p,fl)r<| — l^a 2 ^)^. 



o<p<£ 



ap h' 

Proof. The partial- wave expansion for f (p, 8) is given by 

oo 

f (P'^) = + ^W*) - iJ^CcosW), (6.1) 

e=o 

where a is the radius of the sphere, k = |p, £^ are the Legendre polynomials, and the values 
Sg(K) £ C have modulus one and are equal to 

S,( K ) = = »M±M^, k e R+, (6.2) 

for spherical Bessel functions je, ye of the first and second kind respectively Since the values 
Se(tz) — 1 decay super-exponentially in I and |£^(a?)| < 1 for — 1 < x < 1, the series (16. ip is 
uniformly convergent. It follows the series are also convergent in the norm £ 2 ([0, 7r]; d9sm(9)), 
since the domain [0,7r] is compact. In the consideration of the derivatives of f(p, 8), the 
sums involved will also be uniformly convergent for 8 £ [0,1], since 

Q^^SeOpl) decays 

super- 
exponentially in I and sup_ 1<x<1 |<9™ 2 £^(a;)| increases polynomially for ni,n 2 £ N. 

The identities used in the remainder of the proof related to the special functions je(K), ye(K) 

and £^ can all be found in pp. Some of the facts on the Bessel functions are stated in [I] for 

i i 

the cylindrical Bessel functions J £+ i(k) = (^f) 2 J^(^) and Y £+ i(k) = ( 2 ^) 2 ^(k). 
Part (1): 

Since d z = cos(2~ 1 8)d p — 2p 1 sm(2^ 1 8)de and the factors cos(2 _1 6 ) ), sin(2 _1 6>) and their 
derivatives are bounded, I can simply bound 

/ |p^dpf(p,6>)| 2 for v + w<N, u + w = N, (6.3) 

where u can only be non-zero if v is also non-zero. Since I am considering |p > 1, I can take 
u = v. 

The basic pattern of the analysis will be clear from a treatment of the first derivatives 
dpi (p, 8) and p 1 dgf(jp, 8) along with a few extra comments on extending the bounds to higher 
derivatives. I will focus on the angular derivative first: 

2 00 

p- 1 d e f(p, 8) = -^— 2 + l)(St(K) - l)£j(oos(fl)) sin(0). 
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I have that £,' e (cos(8)) sin(6 l ) = £™(cos(#)) for m — 1, where £™ are the associated Legendre 
functions. The family £™, £ > m is orthogonal for fixed m and have norms given by 

Thus, 

- / |Wp^)| 2 = ^^^+l)(2^+l)|^)-l| 2 (6.5) 



^=1 

For the terms in the sum with £ < [2k\ , I use the crude upper bound 

4 L2«J 4 [2k] 

— i £*(*+i)(2<+ - ii 2 %E^+ 1 )( M + 1 ) 

12™ 4 Z" 2 , o 48yra 4 

where the second inequality uses a Riemann approximation. 

For the sum of terms with £ > [2k\, I require a study of equation ( 16. 2p . From [1], I have 
the asymptotic relations 



„((^ + 2- 1 )sech(a)) = e - (i + 

2(£+ i)(sech(a)tanh(a)) 2 v t J 



a > 0, £ > 1 



+ 2- 1 ) sech(a)) = — (i +0 (^L)). ( 6 .6) 

(£+ |)(sech(a)tanh(a)) 5 V * y 

For £ > \2n\ and fc > 1, 

' < 2e (2<?+l)(tanh(a)-a) < 2 e ~{l+\) 



\S t ( K )-l\ = 2 









for (£ + 2 _1 )sech(a) = k. In the first inequality above, I introduced an extra factor of 2 to cover 
the errors from ( 16. 61) . and the second uses tanh(a) — a < — | for cosh(a) > 2. When k 1 
and £ 3> k = (£ + |)sech(a), then H6. 61) becomes 

j £ ( K )^2-^-|(2£ + l)-5(-^_) £+ ^ and % («) w 2a«"3(2^+ (6.7) 
Hence, for £ > |_2kJ 2 , I have a sharper upper bound than the one above that is given by 

Continuing the analysis for the sum of terms with I > [2k\ , then with the estimates above 



4 oo 

Tra 



H 2 k 



£{£ + l){2£ + 1)\S £ {k) - l\ 



4 



27ra 4 m n , ..wr./, , ,\ -w+i") . 27ra 4 , <WnJ , <w en x «+2 



<^ e ^ + D(« + i).^*' + ^ E ^ + D(« + i)( 5 ^ T ) 

£>|_2/cJ+l ^>L2kJ 2 +1 
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However, the first term will decay exponentially for large k, and the second term will decay 
super-exponentially. 

For the higher-order angular derivatives, (p -1 ^)^ (p, 9) for j = 2, 3, I can follow the same 
routine except I will need to use a few extra identities for Legendre polynomials at the outset. 
For the second derivative, I use that d]j of £^( cos(0)), £ > 2 can be expressed as 

<9 e 2 £ £ (cos(#)) = sin(0)£]'(cos(0)) = 2~\£ + l)££t( cos(0)) - 2~ 1 £ 2 , ( cos(0)) , (6.8) 

where £™' is the derivative of the Legendre polynomial £™. In the above, I have used the 
definition of 2} t ' and the identity 



(1 - x 2 )5£f'(x) = 2-\£ + !)(£- m+ l)£f _1 (x) - 2~ x £f +1 (x), m + 1 > t. 



(6.9) 



By (16. 8p and the inequality ($^ =1 x j) 2 — 3 S"/=i x "j 

p- l d e )H{^6) 

37ra' 



< 



3yra 6 



/*7T v 

/ desm(e)\^2£(£+l)(2£ + l)(S e (K) - l)£*(cos(0)) 



+ 



d0sin(0)| J2(2£+l)(Si(K) - l)£ 2 (cos(0)) 



1=2 



+ 



3W> 

8/1 4 k 6 



(^(/c) - 1) | 2 (6.10) 



where the last term of the right side is the £ — 1 term from the partial-wave expansion. To 
simplify the expressions in (I6.10p . I can then use that the £^(cos(0)) are orthogonal for fixed 
m and have normalization (16. 4ft . For the third term on the right-side of (16.101) . 



3na 6 



d6sm(0)\ J2(2£+l)(S e (K) - l)£ 2 (cos(0)) 



e=2 

3vra 6 

4h^ 



J2{2£+l){£ + 2){£+l)£{£-l)\St(K) - l| 



o 6 L 2k J o 6 



1=2 



2h A K? 



J2 (^ + 2) 5 |S',(«:)-1| 



1= |2kJ+1 



I can apply the same analysis as for the first order derivative by giving a separate treatment 
of the terms with £ < [2k\, [2k\ < £ < [2k\ 2 , and [2k\ 2 < £. It is crucial the degree of the 
polynomial in £ is strictly smaller than the power of k^ 1 appearing in the prefactor. 

The first term in (I6.10p is handled similarly. The third derivative requires two applications 
of (16. 9 p to write (p~ 1 de) 2 £e( cos(0)) in terms of £"(cos(0)) for < u < 3 and is otherwise 
similar. 

Now I bound the first-order radial derivative d p f (p, 0). Differentiating term by term 



d p f(p,9) = ^J2(2£ + l)d K { 



S,(k) - 1 



£=0 



)£^(cos(0)). 



(6.11) 
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By the orthogonality of the functions ££(cos(#)) and the inequality (x + y) 2 < 2{x 2 + y 2 ), 

/ mp^i 2 = ^iE(^ +i )i --(^w - 1 ) +^mi 2 

^ S B 2 ^ + 1) W«) - if + S B 2 ^ + l)|^(-)| 2 . (6.12) 
£=0 £=0 

The first term on the right side is ^^ r cr tot (p), which is bounded over the domain p > 1, since 
°tot(p) is bounded. For the second term on the right of ( 16.1 2p . I can obtain an expression for 
d K S t {K) as 

-2m- 2 



(6.13) 



where h\ (ft) = j^(k) + ye(n)i and I have used the Wronskian identity W[je(K), ye(n)~\ = k 2 . 
Thus, 

9 4 00 o 4 00 —4 

^ E (« +1)MW r=^ E (« + i)^, ,6,4) 

where m v (K) = (y 2 (K,) + j 2 (/t)) 2 . Again I give a separate analysis for the terms £ > [2k\ and 
£ < \_2k\ . For the sum £ < [2k\ , we can use m^n) > k^ 1 for all £ e N and n G M + 

o 4 L 2k J _4 o 4 L 2ft J 00 4 

The sum of terms with £ > |_2kJ can be treated with the estimates (16.61) and (16. 7p . which I 
used for the angular derivative dgf (p, 9): 



<^ l g (2 < +1)e -^) + l^! f (2< + l)(-^-) 4 ', 

1>\_2k\ £>[2k] 2 

where I have doubled the coefficients to cover the error for the asymptotics. The right side is 
exponentially small for large k. 

For the higher derivatives 9^F(p, 9) with w — 2,3, the coefficients in the partial-wave ex- 
pansion include a^(^^i), and it is necessary to find convenient expressions for the derivatives 
d™S(K). Using the relation 

d K hf\K) = -hf\K) - hfl^K) (6.15) 
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and (16.131) . then 

The absolute value is less than 

I 1 2 

Hence, I have the upper-bound for the integral of |(<9pf)(p, 9)\ given by 
|<9pf(P,< 

<^^ (M+ l)( K -«| St ( K )-l| 2 + K - 8 (£ + 3 )V 1 (-)+-- 6! ^). (6-17) 



n 



where I used (x\ + %2 + x^) 2 < 3[x\ + x\ + As usual, the terms in the sum can be split 
into I < [2kJ, \2k\ <£< [2k\ 2 , and £ > [2k\ 2 . For the £ < [2k\ terms, I use 

sup ;r^— < oo. (6.18) 

k>1,0<1<[2k} m A K ) 



This can be seen by using the polynomial expression 

t 



k? ^— ' n!(£-n)! 2 



n=0 



Hence, I can again use m^ft) > for the first \2k\ terms in (16. 17ft . The terms with £ > [2k\ 
can be treated using the approximations (I6.6P and (16. 7p . 

The third-order derivative fit? can be bounded by the same techniques as for second-order. 
I use (I6.15P to reduce d^S^n) to an expression in terms of Hankel functions, and then bound 
\d^Si(K)\ using the m^ys. An inequality of the form (I6.18P will be required for the ratio 

rai +2 (ft) 
mt(K) 

Finally, the mixed derivatives v ^ 0, w ^ in (16. 3p require no new observations from those 
above. 

Part (2): 

I now consider the case k = |p < 1, although my concern is for the regime k < 1, since 
finite values of k are covered by the analysis in Part (1). The factors p _1 appearing in (16. 3 j) 
grow as p — > 0, but this turns out only to prevent a constant upper bound for (16. 3p in the cases 
when u = 3. 

The expression ( 16 .31) for v — 1 and w = is equal to 



r. 2N+2 00 

/ |p-^f(p,^| 2 = ^^^^+l)(2£+l 



SAk) - 1 



K 



]7ia 2N + 2 

3h 2N 



k 2 ( 2 - n )(1 + 0(k 2 )), (6.19) 
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where the first equality is (16. 5ft . and the approximation is for k<C 1 and follows from the power 
expansions 

k 1 / 2~ 1 /t 2 (2- 1 k 2 ) 2 \ 

3 * K ' ~ (2£ + l)\\\ l\(2£ + 3) + 2!(2£ + 3)(2£ + 5) J' 
21-1)!!/ 2- 1 /? 2 (2- 1 k 2 ) 2 \ 

Vi{K) = -^^{ 1 -W^) + 2l(l-2i)(3-2i)----)- (6 - 20) 

With simpler arguments than for Part (1), then (I6.19P has a constant upper bound for all 
< p < - when N = 1,2 and is bounded by a constant multiple of p _1 when N = 3. For (16. 3 j) 
when v = 2,3 and w; = 0, I again use the identities (16. 4p and (16. 9 j) as in Part (1) along with 
the approximation (I6.20p . When «; 7^ 0, then I write the derivatives i9 / ^( ' 5 ' f( ^' > ~ 1 ) in terms of 
Hankel functions as in Part (1) and use ( 16.20p . 

□ 

Part (2) of the lemma below is a technical point, which I use in the proof of the convergence 
of the of the diffusion constant in Thm. 12.11 

Lemma 6.2. Let b, k E R 3 . 

1. There is an asymptotic expansion 

3 1 3 

j=l i,j=l 

for some operators C^\j), C*jP(i,j) = C?\j,i), and that are relatively bounded to 

Cq. Moreover, the constants for the relative bounds are uniform for k,b in any bounded 
region. 

i III II IJU II II (/.(.(J IIAjll/ll 1111:11 II III/ I. T — A. I ll / / 

such that 

\\GyWc* < C n \\y\\ c n+i. 



2. For all b,k in a bounded region and any G = C~ (j), or there is aC n > 



Proof. 
Part (1): 

By Part (5) of Lem. I5.3[ the graph norm of Cq is equivalent to the weighted norm || • ||f^. It 
is hence sufficient to show the operator valued derivatives of £g up to third order are relatively 
bounded to multiplication by \p\ with uniform constants in a bounded region of k. The generator 

splits into three terms 

which operate as in (14.91) . These terms will be treated in Parts (i), (ii), (iii) below for 



(i). h h (ii). J h (iii). S % . 

These operators have domains including 71 by the proof of Part (6) of Prop. 15.31 An implication 
of the analysis here is that the operator derivatives can also be regarded as maps (Vg £g) : 
Ti -> L^M 3 , C 3JV ). 
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(i). The operator £g acts as multiplication with function 

The map V? £^ : 7i — > L^M 3 , C 3JV ) acts as multiplication by the vector- valued function 

(Vf £ % )® = I(?)"((-l)"(V«"£)(jf- + (V*"£)(?+ 2- 1 M)) (6.21) 

Hence, it is sufficient to bound |(V® £){p)\. The vectors (V® £)(p) can be expressed as 




where i^^v is the iVth Hermite polynomial. The second equality is the same change of integration 
variables as in the proof of Lem. 15.11 

Using the above integral expression, |V® £(p)\ is bounded by 




where I have made a change of integration p re \ — > q = y m(ntP rel + MP) ^ or ^ e secon d 
inequality. 

By the triangle inequality, the above is bounded by 

which gives a linear bound in \p\ for (I6.2ip . It follows that V?™££ is relatively bounded to the 
multiplication operator 

(ii). By the form $00]) for J n : Ti ->■ L\R 3 ), 

Vf J % = 4-^ 2 / dqr ? wf\ wf\p) = {Vfj S )(p+ q,p), (6.22) 

where r^-is the shift operator {T^y)(p) = y{p — q) for y G L 1 (R 3 ), and acts as multiplication 
as above. The norm of Vj? operating on an integrable function / is smaller than 

\\Vfjtvh<[ dp\y(p)\([ dq\(Vfjz)(p + q,p)\). 
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My aim is to show J R3 dg*|(Vf J%){p + q,P)\ is linearly bounded by \p\. 

I dq\(Vfj R )(p + q,p)\< V ^f dq^rj dvr(v + 2^^-q + 



fimh z 



e-^ fe ii ? f(Q_,0_)f(Q + ,0 + ) 



(6.23) 



I will take the gradient to have components in a basis with directions q and two orthogonal 
vectors. By the product rule, 

e -^ fc ii 9 -f(Q_,e_)f(Q + ,e + ) 



=e 8m 2 ii? 



E 



AH 



N X \N 2 \N 3 \\\ 4M 2 / ^ 

jVl+jV 2 +Ar 3 =iV 

x (-i)*(vff)(Q_,e_) (vf 3 f)(Q + ,e + ). 

Using the identity (14.111) and the inequality 2\xy\ < x 2 + y 2 , 



(3mh 2 
AM 2 



t 1 171 _, ill v r _ 



<2 : e 8M a "'lis 



E 



f(Q_,e_)f(Q +J e + )]| 

AT! / l/3mh 2 \ Nl 
A^TiVaT/Vs! VV 4M 2 J 



3 



iVi 



/3m/?. 2 -* 
4M 2 " 



x a 



,N 3 -N 2 . 



■(£.) I Vf 2 f (Q_, 0,) I 2 + a^-r(£,) |vf "f (Q +) 6+) | 2 ) . 



In applying the inequality 2|xt/| < x 2 + y 2 , I have multiplied and divided by a factor of a^ N2 
to keep the units consistent. Returning to the integral (16.231) . a single term from the sum over 
N x + N 2 + N z = N can be bounded by 



7/ N\ 



2 N]\N 2 \N Z \ V V 4M 2 



f3mh 2 \Ni 



m 



sup e 37 -Wt))-^ / dq— / dv 



1 



x a 



,N 3 -N 2 



(,-_)|vf 2 f(Q_,e_)| 2 + 



n 2 -n 3 



(i- + )|vf 3 f(Q + ,9 + )| 2 ), (6.24) 



where I have pulled out the supremum for the part of the integrand depending on kn q . 

I will focus on the term r(z_) | VS^f (Q-, 0-) | 2 . It is convenient to switch the integration 
from J R3 d^i / Ra dv to ^ J R3 dp re i J n where 

777< 777/ — * 

Prci = —z- - -rj{p- 2~ 1 hk), so that |p re i| = Q , 
m M 

and the integration over Q reduces to f n = 2tt L , dO_ sin(0_), since the integrand (I6.24[) 
does not depend on the azimuthal angle 0_. The gradient V^f can be written 



Vjfd^i.e.) 



cos(2- 1 e_)(a p f)(|p rel |,6_) -2|p rcl |- 1 sin(2- 1 e_)(^f)(|p rel |,0_) 
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for w = —v— jf(p± q — 2 1 hk±g) and the specific dependence of w on p mh 0_, <j)_ is not 
important here, since only the norms of the gradients appear in (16.241) . More generally, 

vf f(liue-) = (^)"(| f )"'(«m(l«*l.e-), 

where d z is defined as in Lem. 16.11 Notice d z is the derivative with respect to the variable 
z = pcos(2 -1 ^) while holding the variable y = psin(2 -1 ^) fixed. By changing integration 
variables as suggested above and applying Lem. I6.1[ I have the first equality and inequality 
below 

Co 2 [dm* ,2N 2 f w , „ , / m _ m,_ r . N / , H ^2, 

Ca 2 ml ,am^2N 2 f m -i fc A| / 

( w ) y R3 ^K'-M^- 2 «0|r( 



Af,3 



m 4 



w 



where G f2 is a unit vector. The left term on the right side of ( I6.25P can be bounded by the 
triangle inequality and computations with the Gaussian r(w) as in the proof of Lem. 15.11 

The integrand for f n in the right-most term of (16.251) is maximized when the vector p 9 points 
in the opposite direction as p — 2~ 1 hk, so 

/ / dppr(—p9 + —(p-2- 1 hk))<4n dppr(—p-—\p-2- 1 hk\) 

JnJm+ ™* M Jr+ m * M 

I Tfl Tfl ~* 

<4tt / rfp|p|r(— p-— \p-2- l hk\) 

where the second inequality is an extension of the integration domain, and the third inequality 
follows by the triangle inequality after computing the Gaussian integral. 
For constants Cn 1: n 2 ,n 3 > defined as 



77 JV! / f3mh 2 \ N w 12 



r\ iV! / pmh'y^/ _i 2 \ 



m*V v m y a 2 m 3 / vv 7r/3 y M 2M 1 Iy V v 2M ; v 2M 
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then (16.231) is bounded by 

CnxM.Ns < oo. (6.26) 

Nl+N 2 +N 3 =N 

The bound (I6.26j) is linear in \p\. It follows the partial derivatives of [ty]z up to third-order 
are uniformly relatively bounded to \p\ in a bounded region of k G M 3 . 

(iii). The map : 71 — > L 1 (R 3 ) acts as multiplication by 

h n (p) = H(p- T x hk) - H(p + 2- l hk). 

Thus, Vf /i£ is a multiplication operator with function 

(^) iY ((-l) 7V (V^^)(p-2- 1 M) - (V 0Ar //)(p + 2- 1 ^)) 

= -i^4p + i(^) 7V ((-l) JV (V^i/ / )(p-2- 1 ^) - (V® N H f )(p + 2-ihk)), 

since = fjg-p 2 + Hf(p) where Hf(p) is defined in (12. 6p . By a change of variables, I can 

write 



Hf(p) 



The derivative has the form 



mz 



/ dp rcl r ( — p ml + —p) Re [f ( \p TC y \ , 



Hf)(p) 



2nrjh m ( j (3m\ 



mz 



N 



dp 



rcl 



m _ m . 



x r(—p vel + ^-p) Re[f (|p ro i|, 0)] . 



By the same change of variables as for (i), the term |(V® Hf){p)\ is bounded by 

1 



2-KT]h I m i / U3m\ N 
~ A W At*) 



mz 



dp 



rcl 



N 



< 



< 



2jrh 2 r] / U3m\ N 
1 Jp) 



m _ m _, 
m* M 



dq\Re[f(\(j)^q- ^,0)]||MI5I) 



r (^~^ rel + Imp) l Re [ f °)] 



2'ircK 'q / /3m\ N 
1 MV 



#(£ + (j)*l*l + ^0)Ml* 



e 2 



e 2 1 



(27T) 



(2tt; 



3 " 

2 



(6.27) 



where the last inequality uses |f(p, 0) | < c{a 1 h + p) for some c > 0, which is shown in the 
proof of Lem. I0.lt However, the functions p3jv(|<z| 



and \q\\Sj N {\q 



-y 2 



T are in L 

(2tt)2 



as used in Part (i), so (16.271) gives a linear upper bound for (V® Hf){p) \ in terms of \p\. The 
operator V? /ig is therefore relatively bounded to multiplication by \p\. 

Part (2): 
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By Part (4) of Prop. 15.31 the norms || ■ ||£« and || ■ are equivalent, and it is sufficient 
to prove HGyHa^n < C||y|| OT n+i for G = V? £g, N = 1, 2, 3 and k in a bounded region. This 

will be convenient, since M. n commutes with the parts of V- £g corresponding to the loss and 
Hamiltonian parts of £g (as in (i) and (iii) of Part (1)). By the triangle inequality, 

\\Gy\\ mn < \\Gy\U + \\GTTy\U + \\ [Wl n , G]y\\ 1 
< c\\y\\ m n +1 + \\ [mt n ,G]y\\ 1 , 

where c > exists by the results of Part (1) and the equivalence of || • ||£ and || ■ [| sxrr- 

Since DJl n commutes with multiplication operators, [271", Vf "jCj] = [371", Vf"^]. By (E22D, 
the commutation [9JT n , V? can be written as 

m n ,Vfji=4-w[ dqr,zf\ 

where is multiplication by 

zf\p) = (\P+ ?T - m<yfj S W+ Q,P)- 

Again I have 

||[9JT\Vf Jdv\\i< [ dp\y(P)\ [ dq\zf\p)\. 

The value | z~* (p) | is smaller than a product of | \p + q\ n — \p\ n \ and | V? J%{p + q, p) | , and by 
(ii) of Part (1), 

N!+N 2 +N 3 =N 1 z 6 7fcK + 

x S W I ^(^ 3 -" 2 K--)|vr 2 f(Q-, e_) | 2 + a— 3 r(£+ ) i yf f (Q +> e + ) | 2 ) . 

For factor \\p + q\ n — |p| n | has the brute upper bound 2 n (|p| n + |g| n ). Applying the stan- 
dard change of coordinates from q, v to p TC \, 6 and following the reasoning in (ii) of Part 
(1), then J R3 dq\z^\p)\ can be bounded a degree n + 1 polynomial in \p\. By interpolation 
||[9JT n , V? is bounded by a multiple of ||y||x n+1 - 

□ 



7 The diffusion constant 

The following proposition relates the diffusion constant D = _Dkm + -Dj P s to the second derivative 
for an eigenvalue for C^. Part (3) is another another technical point required for the proof of 
the convergence stated for the diffusion constant in Thm. 12.11 

Proposition 7.1. Let g > be defined as in Prop. 15.31 There is an r > such that for all 
\k\ < r, the following statements hold: 
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1. The generator has an isolated eigenvalue e(k) G R_ + iR with 

e(k) = - 1 -D\k\ 2 + 0(\k\ 3 )- 

2. The spectral projection corresponding to e(k) is non-degenerate and has a Taylor ex- 
pansion 

P % = P + ki(%Pn\**) + Yl hkjididjP^o) + \k\ 3 m f (k), 

i i 

where 9jP^|^ =0 , didjPz\z =0 and 9i'(k) are uniformly bounded with respect to the operator 
norm for linear maps on L 1 (IR 3 ). 

3. For n > 1, there is a C n > such that for any G = diP$\$ = Q, didjP#\ff =0 , or 9t'(k) 

\\Gf\\ C n < CjfWjq-L 

4- The spectrum of satisfies — {e(k)} C (— oo, — 2~ 1 $j] + iR. 

5. The semigroup e tC k with generator C-. = + § (I — Pg) (having domain D(£y = D(£^) ) 
satisfies 

\\e tC 'n\\ < C 

for some C > and all t G R + . Also, 1 1 e tc k (I - P^) 1 1 < Ce~* * . 

Proof. 
Part (1): 

By the second-order expansion in (16. 2ft with operator derives relatively bounded to Co, 
perturbation theory [27] guarantees there is an isolated eigenvalue e(k) near zero with expansion 

e{k) = k-3L + 2- 1 k® 2 • (b- c) + 0{\k\- 3 ), 

where a G C 3 , b, c G C 3 £§> C 3 are of the form 

a, = {1\C!PU)M, b M ' = (l|4 2) (^i)koo), c id = (l|4 1) (i)S'4 1) (i)|x/ 00 ), 

where 5 is the reduced resolvent of Cq. Note that (l^l^oo) = J K3 dp{Av 00 ){p) for an operator 
A on L 1 (R 3 ). I claim a = 0, bjj = —di^D^, and c^.j = <5jjDkin- 

I begin by finding expressions for the first two derivatives at k = 0. The derivatives are 
given by 

(V^U = -ih~\V® N h f )\ p=0 + (V^)|, =0 - (V^^)|, =Q! N = 1, 2, 

where the terms on the right are discussed presently. As a result of Lem. 16.21 the operators 
above map T x into L^R 3 , C 3 *). 

The operator — (V/v)|r=o ac ts as a vector of multiplication operators with multiplication 
function h(VH)(p). The terms (V® h r )\ F=0 and (VSf)\ r = o are both zero due to symmetry 
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around r — 0. To write (V® 2 £?)|f=o, (Vj7r)|r=o, and (V® 2 J7f-)|f=o, I will make use of the func- 
tions Tg{pi,p 2 ) in the form ( 12. 3ft . The operator (V® ^l-.^ acts as a vector of multiplication 
operators with multiplication functions 

{V*S?m\^ = 4T x tf [ dq((V 1 + V 2 r 2 T ? Mp), (7.1) 

Jr 3 

where Vi and V2 are the gradients with respect to the first and second arguments of T^(pi,p 2 ). 
The operators Vj7^|^=o and V® JTr\r=o have the form 



(VJf)| F=0 = 2-'h / dqr^, (7.2) 
Jr 3 

I -1*2 f j^„.,(2) 



(V* » M = 4-^ 2 / dgr^, (7.3) 

Jr 3 

where acts as a shift {r^y){p) = y(p — q), y e L 1 (M 3 ) and and tu^' are vectors of 
multiplication operators with 

4 1} ($ := -(( V i - V 2 )T g -)(p,^) and wf{p) := ((Vi - V 2 )® 2 T g -)(p,p). 

Now I am almost ready to evaluate a, b, and c. As a preliminary observation, there are 
rotational symmetries 

u 00 (Rp) = u 00 (p), H(Rp) = H(p), T Rq iRp u Rp 2 ) = T q ip u p 2 ), R e SO s . (7.4) 

It follows that a must be zero. Moreover, the tensors b and c must be invariant under the 
operation of R £g> R for all rotations R, and therefore b, c are constant multiples of the identity 
tensor. 

With ( 17. ip and ( 17.31) . the value for b can be written 

b = <l|(V^U- V ^H,=o)^> 

= -h 2 dp dqu^ipjiVi ® s V 2 T f )(p,p), 
Jr 3 Jr 3 

where £g> s is the symmetrized tensor product. By the rotational symmetry mentioned above, I 
can write 

[(IKV^JfI^-V^I^)^)] =-6^3-^ [ dp I dqu^AT^p), 

Jr 3 Jr 3 

where A = Ylj=i 2 3 ^ij^.j; and 9jj is the derivative for ith component in the ej direction for 
some orthonormal basis (ei,e 2 ,es). The definition of A is invariant of the basis used. The 
expression (AT^){p,p) is equal to 



r 3 

(AT»(£p) = / dvY" \djL^(p)\' 



m 1 

--7] 



/ dvr(v + 2 — q+ —p M ) 

Jin), m * M 



™l \q\ j %A 
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where p re[ = ^v-jfP±^+2 l q and 2 tan(2 1 Q) — |g| | ^ u— ^p± g | . In the above, the basis 

vectors (ei, e 2 , e 3 ) for A have been chosen to depend on q, v,p as e 1 — e 2 = -=-h-, and 

with e 3 as one of the two normalized vectors orthogonal to e\ and e 2 . The derivative d^L^{p) 
with respect to the third direction e 3 is equal to zero. 

By using equality (17. 5p and by shifting to center-of-mass coordinates for the integration 
variables 



dp / dq (p) ( ATq-) (p, p) 



g 2(m+M)rcm P g 2 m» ^rel 

'cm 3" / *^Prel |.Prel | 3" 

(27r(m + M)/?- 1 ) 5 (2nm^- 1 ) 5 

?2 ^ ™2 



X 



= / rfprel Prel - ~T ( TTTJ (iPrelT + 3 — — ) a tot (|p re l|) + T7^ z (|p re i|) ) 

J R3 (2Tim 1r j3- l Y V4M P ' 

I have shown 

= [<i|( - j^MU + (v^)U - (v« 2 ^)|^ )^>] . . 

4r]h 2 /2m M | /" 3 2/ /? , 2 , 2m* 1 » // 2m *a \A 



which is equal to — 5i 5j Z?j ps . 

Finally, to find an expression for c, I again use the rotation symmetry to get 



cm = 3- 1 ^E(ii4 1) h^4 1) hio- 



n=l 



For the image of the spectral projection I — P corresponding to the values S(£ ) — {0}, the 
reduced resolvent is equal to S = — / °° dt e tCo . The image of I — P is the set of elements 

in L 1 (IR 3 ) with integral zero. My earlier observation that a = means (II^q 1 ^)!^) = 0, 
n — 1, 2, 3 so £^(11)1/00 is in the image of I — P. Hence, the elements of c can be expressed as 

6 n=l Jo 

= -^y) / ^(<l|0ne* £o o n |^ 00 ) + <l|o n e* £o m n |z/ 00 )) 

= ^fJ2 ^(llone^o^z/oo) (7.6) 
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where t>„ and ro n are multiplication operators acting on L 1 (IR 3 ) whose corresponding functions 
(denoted by the same symbol) are specified below. The functions t> n , to n are the nth components 
of t), to : M 3 — > M 3 for D defined as in (12.141) and 

1 h r - 



to[p)\= 7^7 77 / dqwyip- q)Voc(p- q) -*>(p)- 

Vco{P) 2 J R 3 

The second inequality in (17. 6p follows from the relation ^l|£Q^(n) = — i(l|t) n and the definition 
of ro. The function tv(p) turns out to be zero by the detailed balance symmetry, since 



K> if) =77 / dq I "™ P Q wf{p- q) - w^lip)) 



' I dq^j Im((Vi7q)(p- Q,P~ 0) ~ (Vi7_ g )(p,p) 

Ju 3 v 2 (n\ V 



2 



V$o{p) 



=0. (7.7) 

where Vi is the gradient of the first argument p\ of j q (pi,P2)- The first equality in (17. 7h is by 
the equality o(p) = | J R3 dqw^\p). The second equality in (17. 7p follows from 



(1) (p) = -((Vx - V 2 )T f )(p,p) = -2iIm[(V 1 T f )(p ! p)] = -2i ^ Im[(Vi7 g )(p,P 



where the second equality is T^(pi,p 2 ) = T^(p 2 ,pi) (since Tq- is the kernel for a positive opera- 
tor), and the third is the detailed balance form (I4.12p . The last equality in (17. 7p holds because 

(Vi7,)(p - q,P~ q) = (Vi7- g )(?,p) by flUB^ 

By working with the adjoint semigroup e ££ «, I can write the elements of c as 

Ci >i = J rJ2 / dp(e tc h n )(pjv n (p)v^(p) 

dm Voo [o(p t ) • ti(po) 



3 

where p 4 is a Markov process whose probability densities evolve according to the master equa- 
tion (12.71) starting from the stationary state z/qo. 

Part (2): 

This follows from the relatively bounded perturbation expansion from Part (1) of Lem. 16.21 
Part (3): 

I will build up control for the Taylor expansion of through control of the resolvent c }_ 

k 

for z in a certain region of z bounded away form the origin. More precisely, I will use that the 
projection can be written as a complex integral 
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for a counter-clockwise curve 7 along the circle of radius 2 _1 g around the origin. 

By Part (1), I can pick r > small enough so |e(fc)| < 4 _1 g for all \k\ < r and £(£ju) 
contains no other values in the ball of radius |g around the origin. Let us also pick r small 
enough so sup|£| <r ||(£g — £o)y\\i < $\\y\\co f° r & U V G L 1 and some small 5 > 0. The resolvent 
expansion 

Ce — z Lq-z^ \ Cq-zJ 

ft 71=0 

makes sense for \k\ < r and |g < \z\ < |g, since there is the following operator norm bound: 

1 



1 



£0 - z 



=5 



£0 — z Co — z 
<S(2 + -C S + - + -C, ) < 2 _1 . 

V 2 5 g g 8y - 



The second inequality above is for some C| > by Part (4) of Prop. 15.31 and the last inequality 
is for small enough 5. 

The first two derivatives for the resolvent around k = are 



0, 



1 



*~"V Z 



\v=0 



v=0 



' ■£?»(.•)- ' 



£0 ~~ 2 Co — z 

1 -cfHij) 1 



Co - z 



where £o(i) and £q (i, j) are defined as in Lem. 16. 2} and the sum is counted twice with % = j. 
The third-order error EjP(z) for expanding c }_ z around k = has the form 



*f(z) : = -A -A k.V(-A— 

k *--k~ z £0 - z Cj-z 



v=0 Cit — z 



v=0 



1 



-R 



(3) 



1 



1 



-i2 



(2) 



(1) 



Co - z k C - z Co - z k C - z k Co - z 



-R 



(i) 



R 



(2) 



Co - z k Co - z k Co - z 



-R 



(2) 



R 



(2) 



Co - z k Co - z k Co - z C 



1 00 

— D-i) B ((^-A,3 

' n=3 



1 \" 



Co - z) 



where R^f 1 is the error for the (n — l)th order Taylor expansion of £g around zero: 



R 



(n) 



c s - £0 - £ ^ ■ (-J— ) 

Cij — z' 



v=0 
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The technique for bounding the second derivative and the error El (z) will be clear from 
the analysis of the first derivatives di c }_ z | g . For \z\ = 2~ 1 $j, I have the following string of 
inequalities. 

II (di£^- z \*Jy\\c« < 4 (! + + c i )ll4 1) « z ^?/ll^ 

<4(l + I)(l + Cj )C B ||_!_ y || 

< 16(1 + -) 2 (l + c £ ) 2 C n ||y|U-i. 
4 

The first and last inequalities hold for some ce by the bounds on II 7^— II and Ht^-H from Part 
(4) of Prop. 15.31 The second equality above is by Part (2) of Lem. 16.21 Since the integral (17. 8ft 
is over a compact curve and all the estimates are uniform, I obtain the necessary bound on the 
value of ||(<9jPtf|tf =0 )2/|U«. 

Part (4): 

For any ball B r of radius r around the origin, there are constants a, b > such that 
sup Il^0lli<a|l^o2/||i + 6||2/||i, yeL\R 3 ). 

\k\<r 3=1,2,3 

Consequently, for < r, the operators C^ — Co are relatively bounded to Co with uniform 
constants given by 



\\(C^-C )y\\ 1 < f dv\\kVC v% y\\ x 
Jo 

< 3ar||£ 2/||i + 36r||y||i. 



I can pick r to make the coefficients 3ar and 3br arbitrarily small. In particular, I can pick r 
so that 

S(£ g ) c E(£ ) + B s , 

and the non-degenerate eigenvalue e(k) is the only spectral value within a radius 5 of 0. It 
follows the rest of the spectrum has real part less than — q + 5, and I can take 5 < 2 _1 £j. 

Part (5): 

By similar reasoning as for Part (3) of Prop. 15.31 there is c > such that for all t G M + and 
\k\ < r for r > from Part (4), 

||e*%(l-P s ) || <£*"**. (7.9) 
The difference between Cr and C'^ is bounded, and thus C'^ is a closed operator with domain 

i K k k 

D(C£). The semigroup e tc * can be expressed as 

e tc '% = e f l e ^(l - Pj) + e u ®P n . 
By the triangle inequality and (I7.9p . 

||e^|| < e^||e* £ s(I-Pg)|| +e te ^||Pj|| <c+l. 

□ 
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8 Proof of main result 



With the results from previous sections, the proof of the convergence for the rescaled position 
distribution in Thm. [2~T1 has the pattern found in [3]. The convergence of the moments requires 
more effort due to the unbounded operators in this case. 

Proof of Thm. \2.1\ Parts (i) and (ii) below treat respectively the distributional convergence of 
the measure t^fj, t (t^r) as t — > oo and the convergence of its second moments. The notation 
dj, j = 1, 2, 3 and V will refer to partial derivatives and the gradient, respectively, of the fiber 
variable k E M 3 . 

3 1 

(i). To prove the distributional convergence of the probability densities p t (t2r) to a Gaussian, 
I will prove the pointwise convergence of the characteristic functions 

<p t (g) = [ drti^r)j™ 
Jr 3 

to e~"?^\ The characteristic function for p t can be written 

ik(k)= f dxp t (x)^ = Tr[p t e^] 

JR 3 

= / dpp t (p,p + hk) = / dp[p t }^(p) 
Jr 3 Jr 3 

= [ dp(e^[p] n )(d = (1| e^[p] % ), (8.1) 
Jr 3 

where I have used that p t e lkX is a trace class with integral vernal K,(pi,p 2 ) = Pt(Pi,P2 + Kk) for 
the third equality, a change of integration variable p + 2~ l hk — > p for the fourth, and the fiber 
decomposition for the fifth. The last expression is the evaluation of the vector e tc% [p\^ € L X (R 3 ) 
for the linear functional determined by integration against the constant function 1^3 G L°°(IR 3 ). 
Fix s G M 3 and let k = t~*s. By the change of variable t^r—> r, 

ip t (s) = [ dxptixY 1 '^ 3 = fl t (k) = (1| e^[p]f), (8.2) 
Jr 3 

where the last equality is from ( 18. ip . 

Now I study the large t asymptotics for e tc k[p]^ e L X (]R 3 ). By Prop. 17. 1[ for small enough k 
there is an isolated non-degenerate eigenvalue e(k) with projection such that the remainder 
of the spectral values have real part falling below — 2~ 1 q. The operator e tc k can be written as 

where the operator norm of e tc z(l — P^) decays with the exponential rate e~i l by Part (5) of 
Prop. EH Thus, by Parts (1), (2), and (5) of Prop. ITTTl 
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where the error is with respect to the || • ||i norm. By Prop. IB.ll the condition Xjp, j = 1, 2, 3 
(and equivalently pXj) is trace class implies [p]% has a bounded derivative in k in the L 1 (M 3 ) 
norm. In particular, [p]^ = [p)o+ 0(t~z), and it follows 



e t£ *[p]g = e-T s " I / 0O +0(H), (8.4) 
since (l|[p]o) = J R dp[p] (p) = 1. Plugging ( 18. 4p in to (18. 2ft . then I conclude 

^(5) = + 0(H). 

3 1 

This proves the convergence of t^p t ^^) t° a normal distribution with variance D. 

(ii). To show the second moments converge, I can rewrite the variance in terms of e s [p]j£ as 

t~ lr Ti\X i Xjp t ]=t~ 1 J dxXiXjp t (x) 

= -r l d i d j (i\e^[ P ] % )\^ 

+ rV^-^ eWi)(I - Pg)[/%>| & - =0 , (8.5) 

where the second equality follows from differentiating (18. ip . and the third is (18. 3p . However, 
the differentiability of [p]^, e(k), and P^ gives 

t-%d,(l\ e^Pip}^ = S^D+Oir 1 )., 

since Ve(fc) = 0. 

The remainder of the proof is concerned with bounding the last term in (18. 5p . Define 
CL = £ % + 2-^(1 - Pj) as in Part (5) of Prop. O Then similarly to (El, e* £ 's = e* e ^Pg + 
e'( £ E + 2)(I — Pj.). By the triangle inequality 

e^ + i)(I - P^W 5 >|^| = \dA(M (e C * ~ e^P^ph^J 

< II^^Wril^Hi + ll^a^e^P^^I^Hx. 

The rightmost term is finite by the comments above. I claim \\didj (e tC k [p\A |^ =0 ||i < c(l + t 2 ) for 
some c > 0. This would complete the proof, since the last term in ( 18. 5ft has the exponentially 
decaying factor e~*2. 

By the product rule, there are four terms to deal with: 

di^(e"%y |g =0 =(^e t£ s| jE L =0 )|>] + (die C '%\ k - ) (djlp]^) 

+ (^%=o) {dilpUtJ + ^ a (didAphkJ > (8-6) 

where the middle two terms are similar. The formal expressions for die tC k\^ =0 : 7i — > L 1 (M 3 ) 
and didje tC k\j: =0 • T2 — > L X (]R 3 ) are given by (18.81) and (I8.10p . Part of the analysis involves 
showing \\die £ |£ =0 ?/||i an d W^idje £ |£ =0 ?/||i are bounded by the multiples of the norms \\y\\m 
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til 1 

and || y Hon 2 5 respectively. The existence of the first and second derivatives of e * is established 
with the third-order error (18.111) analysis at the end of the proof. 

Now, I bound the terms on the right side of (18. 6p . By Part (5) of Prop. 17.11 for all k G M 3 
in a small enough neighborhood around the origin and all t > 0, there is a C > such that 

||e*^||<C. (8.7) 

In words, e i is a bounded semigroup. The last term in (18. 6p is 0(1), by (18.71) and the 
differentiability of [p]^. To bound (0i e s |j£ = o) (dj[p]j:\^ =0 ) in L 1 (R 3 ), note the derivative of e tC k 
at k = has the form 

tc' 



where 0^1^ = d^lfco ~ 2 'fl ^P^. 
For y G L^R 3 ), 



H^ls^yll! < C /"cftill^l^SIL 
Jo 

<C [\t 1 \\d i C^ =0 e t ^h\\ 1 + 2- 1 Cg f dt^P ^ k=0 
Jo Jo 



<CC / dt 1 \\e t ^y\L + S-'Ccofl / dtille^H 
Jo Jo 

< tC 2 (C + 2- 1 c s)\\y\\ £o < tcC 2 (C + 2- 1 c o0 )|| 2 /|| OT . (8.9) 

The first inequality above uses the bound for e tCo above, and the second is the triangle inequality. 
For the third inequality, I use that there some C > such that ||<9i£g|fc =0 2/||i < Co II 2/ II £ f° r a ^ 
y by Part (1) of Lent. 16.21 and maxj ||0iPr|? =o ||oo = c o < 00 by Part (3) of Prop. 17.11 The fourth 
inequality follows because e tc ° and Co commute and then by ( 18.71) again. The last inequality 
is || • \\c < c ll ' Hot where c exists by the equivalence of the norms || • ||£ and || • | ] sot from Part 
(5) of Prop. 15.31 With the above, 



me% =0 )(dAphk =0 )h < tcC 2 (C + 2- 1 c g)\\{d 3 [phk = o) 



but 



= tcc 2 (c + 2- 1 c o0 )(|| {dMuk=o) 111 + II \p\ (PMsls^) lli) 

3 

< tcC 2 (Co + 2-^00) (2- 1 1| {X 3 ,p}\\ 1 + 2- 2 II {Pn, {X 3 ,p}}\\x). 



n=l 



The second inequality is by Prop. IB.lt an d the right side is finite by the assumptions on p. 
For the (c^e* £ |- =Q ) [p] term, 

Jo 

V fdt 2 / t2 ^ 1 e( i - i2 ^'(0n 1 /:i A : = o)e {t2 - tl)z: «(0n 2 /: , |fe = o)e' l£ », (8.10) 
,,. .\ Jo Jo 



+ 



( ™ )= {(K) 
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where the sum over (ni,nz) counts twice when i = j. The first term on the right side is 
relatively bounded to OJt by the same argument as for ( 18. 8ft . Also by the steps in (18.9ft up to 
the next to last inequality, then a single term from the sum in (18. lOj) can be bounded by 

f dt 2 rtftie^^^JjI^e^^^iCy^e^y 
Jo Jo 

<C 2 (C + 2- 1 c oQ ) [ dt 2 [ * dt^djC'^y^'oyL 
Jo Jo 

<C 2 (C + 2- 1 c g) [ dt 2 [ 2 dt 1 (C \\e t ^y\\ c2 + 2- 1 c fl||e tl S||i) 
Jo Jo 

< 2- 1 cC 3 t 2 (C + 2- 1 c o0 ) 2 || 2 /|| OT2 . 

The second inequality follows by using the triangle inequality with 9jC'A^ =0 = djCj^\^ =0 — 
2 -1 g<9jPg|g =0 and applying respectively Part (2) of Lem. 16.21 and Part (3) of Prop. 17.11 as 
before. 

Since I have used the formal expressions ( 18. 8 j) and (18. lOj) for the derivatives of e & at k = 0, 
I should control the resulting error in using those terms in a second-order Taylor expansion. 
Using a rearranged third-order Duhamel equation, 

Ef :=e tC '% - e tCo - k ■ (Ve"£|*=o) - 2" 1 fc® a • (V®V^|f=o) 

dt x e^nm e ^o 



+ 



f dt 2 dti e^^R^e^^'oR^e^ 

II k k 

Jo Jo 

f dt 2 f 2 dt^-^R^e^^Rfe 1 ^ 

Jo Jo k k 

f dt 2 dt x e^^R^e^-^'oR^e^o 

Jo Jo k k 



+ / dt 3 
'o Jo 



t-2 



dt 2 / ^e^^^^e^-^^o^We^-^OA^a)^^ 

' k k k 



(8.11) 



where R^ is the error the nth order Taylor expansion of C'~ around k = 0: 



4 n) = 4 - C'o - E ^ ■ ( V ^k=o), n > 0. (8.12) 



m=l 



(n) 
fc,2 



Since £' s = £r + 2~ 1 0(I - PA, the error #l n) is equal to P^ - 2- 1 gi?^ where P^ and R 
are the analogously defined errors for £g and P^, respectively. 

Using Lem. 16.21 and Prop. 17.11 along with the same techniques as above for the first two 
derivatives, then it can be shown 1 1 El^?^ 2/ ] 1 1 is bounded by a constant multiple of |/c| 3 ||y||OT 3 for 
all y e T 3 . By Prop. EU 

\MzWw < ll[p]felli+ WPhPiJpMih 

ii,i2,i36{l,2,3} 

<IH| 1+ ]T 2- 3 ||{P il ,{P, 2 ,{P i3) P}}}lli- 

^£{1,2,3} 
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By the triangle inequality and the self-adjointness of p, the term on the right is bounded by a 
sum of terms HP^P^ppJIi and WP^P^P^pWi. However, 

HP^pPjIl < || |P|V%lll = UPaPl^Hll < II I^IPl^Hll < OO. 

The first inequality holds since {P^pP^P^ \ 2 P i2 pP i3 )^ < (p 3 p|P| 4 pp 3 ) 2 by the operator mono- 
tonicity of the square root function. The equality uses that the adjoint operation preserves trace 
norm, and second inequality is by the same reason as for the first. Finally the right side is 
finite by my assumptions on p. 

□ 
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A A Fourier transform relating ^ and the component 
maps 

In this section, I point out a canonical relation between the completely positive map \I/ and 
the maps T^- in the form (12.21) through a Fourier transform. The statement and the "proof of 
Prop. IA.1I are heuristic. I think the basic mechanics for the relations may interest the reader, 
but it is not necessary to clarify the analytic details for this article. 

Note that the covariance properties flA.2j) below is equivalent to saying the maps : 



B\ (L 2 (IR ci )) act as multiplications operator in the momentum representation as in (12. 3p . 
Proposition A.l. 

1. Let \I/ : £>i(L 2 (R d )) be a map of the form 

tf(p) = f dqe^T^e-^, (A.l) 



where the maps are completely positive and for all x,q EM. d satisfy 

T^rgp) = TgT^p) and T^pr s ) = T^p)^, (A.2) 

in which t$ = e l T> p . Then, \I/ is translation covariant and completely positive. Moreover, 
the map is related to the maps through the Fourier transform 
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2. Conversely, if \& is completely positive and maps Tg- are defined by $A.3\) . then the maps 
Tg- are completely positive and satisfy $A.2\) . Also, ^ can be written as ( L4.il) . 

Proof. 
Part (1): 

The complete positivity of \I/ follows trivially from the form flA.ip . and the translation 
covariance of ^ follows from the the left and right covariance properties (IA.2P and the Weyl 
intertwining relation 

eHMVf* = e^e^e^. (A.4) 
Also due to the intertwining relations ( 1A.2j) and (1A.4j) . 

= T f (p). (A.5) 

The above comes through switching the order of integration and the formal identity 5(p — q) = 
i f dre^^-P) 

Part (2): 

First, I show satisfies (1A.2j) . By the Weyl intertwining relation (IA.4j) . 

= e "' fx X(2^/„/ feifT i*<^ rf )^ 

= T/np), (A.6) 

where the second equality is a change of integration x — y — > x, and the third is the definition 
of Tg-. The same argument can be applied to show T^(prg) = T^(p)rg by using the expression 
r^{r x p) in the integral ( 1A.3j) can be replaced by ^{pr*)rg due to the translation covariance 
of . 

Now, I show Tg- maps positive operators to positive operators if \I> is completely positive. 
The argument is based on Bochner's theorem. First notice the translation covariance of \I/ 
implies T^(TZpr y )r* G 23i(L 2 (M d )) is a function olx-y. Let G L 2 (M d ), p G £i(L 2 (M d )) be 
positive, and p : R d — > C be defined as 

Fj,, P {x-y) := ((p\rg^{r^pT^Ti(j)). 

By the complete positivity of it follows F^ iP {x — is the integral kernel for a positive 
operator on L 2 (M. d ) when p is a positive operator. By Bochner's theorem the Fourier transform 
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of F^p is positive-valued, so 

< / dx e'T F^ p (x) = \4> ( dxe l ^r*:^{T S p) 

JR d ^ ^ JR d 

= <e- i ^0|T,(p)e- i ^0>. 

Since <p £ L 2 (M. d ) is arbitrary, it follows that Tq(p) is a positive operator. The argument above 
extends trivially to show that is completely positive by replacing \1>, T^- with \l/ Cg> I ext , 
T^® / ex t operating on the extended Banach space Bi(L 2 (M. d ) <g> H cxt ) for some Hilbert space 

A computation similar to HA.5[) shows that plugging in to the right side of (lA.ip yields 



□ 

B Properties of the fiber maps 

The following elementary proposition gives sufficient conditions for the smoothness of the fiber 
elements [p]^ p £ B 1 (L 2 (R d )) as a function ofkeR d in the L 1 (R d ) -norm in terms of weighted 
tracial semi-norms of p. The proposition also gives bounds for the moments of [p]%. The partial 
derivatives dj, j £ [1, d] in the proposition below refer to the fiber variable k. 

Proposition B.l. Let p £ B 1 (L 2 (R d )), then 

1- \\d n ■ ■ ■ d jn [p)%h < Mi X hr ■ ■ > i X ^P) — >IU, 

3- IIM[p]dli<ill{^{^p}}lli- 

Proof. By the definition of the fiber maps : Bi(L 2 (R d )) — > L 1 (R d ), they satisfy 

dp gift IpUp) = Tr g(P + 2- 1 hk) P ] 

= Tr[g(P) e^ x p^ x ] 
for all g £ L°°(R d ). Differentiating both sides by dj and taking the absolute value, 

dpg(p)dMM < 2- 1 |Tr^( J P)e i ^{X J ,p}e i ^]| 
<2- 1 ||( 7 || 00 ||{X J ,p}|| 1 . 

Supremizing over all ||^||oo = 1, the left side is equal to Higher derivatives work by 

the same argument. In order to bound the moments in Parts (2) and (3), recall e tkX pe lkX = 
e lk { x ' }(p) and {Xj, •} and {Fj, •} commute for all 1 < i,j < d. 

□ 



/ 

Jr° 
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C Some relative bounds 

Lemma C.l. The operators Hf and are relatively bounded to \P\. 

Proof. In the expression (12.51) for L^$, the only unbounded factor for fixed q E M 3 , v E (q)± is 
the scattering amplitude 

= f f I"- - ^+ 2"^, Stan- 1 ( , 2 "'^ , )). 

I to M ^-L9| 

In particular, the unboundedness of f (p, 0) occurs as p — > oo for 9 near zero. However, there 
is a c > such that 

|f(p,0)| < c(a- 1 /i + p). (C.l) 
If I show this, then is relatively bounded to \P\, since 

VmM m M / v m M ' 

Showing (IC.lj) requires a similar form argument as in Lem. 16.11 By the partial-wave expan- 
sion (16. ip and the fact that |P£(cos(0))| < 1, then 

- oo 

i f (p^)i<^-E( 2£ + 1 )l^^p)- 1 l- 

However, for |p 3> 1, this sum is bounded by a constant multiple of p by the analysis in the 
proof of Lem. 16.11 

For the Hamiltonian Hf, 

,2, 



■ / dqr(q) (a~ 1 h + \—q- -r^-p]) 



2-Kch 7] 



m 

^2 



2irch 2 r) ( _ x ra* ,2m,\ m* 

fo h + 2 — (— t) 2 +77 

m* V m up M 



where I have used inequality fIC.ll) . an explicit integration of the Gaussian r(q), and the triangle 
inequality. The above gives a linear bound in \P\ for Hf(P). 

□ 
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